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ABSTRACT 

A discussion  is  given  of  ths  physical  ant!  jaathem&tieal  background  of  the  wave 
phenomena  produced  in  materials  by  periodic  variation  of  the  tsnpers.turs  or  heat 
flux  at  a boundary.  It  is  shown  how  the  corresponding  mathematical  problems  separ- 
ate into  o superposition  of  a static  and  a dynamic  problem,  where  the  former  satisfies 
Laplace"  - equation  and  the  latter  the  wave  equation  with  a complex  propagation  constant 
but  with  simpler  boundary  conditions.  In  the  case  of  the  dynamic  problem,  a rigorous 
transmission  line  description  is  developed  for  the  principal  mode  temperature  waves  travel- 
ling on  cylindrical  rods,  where  the  voltage  and  current-  correspond  respectively  to  the 
principal  mode  temperature  (constant  over  the  cross  section)  and  the  total  heat  flux 
through  a cross  section.  Obstacles  or  discontinuities  in  the  rod  are  then  exactly  repre- 
sented by  lumped  circuits  on  the  transmission  line,  provided  that  communication  between  • 
obstacles  or  sources  takes  place  only  via  the  principal  mode  with  its  considerably  lower, 
although  still  finite,  attenuation  compared  to  the  higher  modes.  It  is  also  necessary 
that  th®  characteristic  Impedances  of  different  transmission  lines  meeting  at  a junction 
satisfy  a particular  condition  in  order  that  the  junction  shows  the  proper  reciprocity 
properties  characteristic  of  a circuit.  This  condition  is  shown  to  be  satisfied  by  the 
characteristic  impedances  resulting  from  the  above  definitions  of  current  and  voltage. 

A general  theorem  restricting  the  va  ue  of  any  thermal  impedance  to  be  resistive  and 
capacitative  ( i.e.  restricted  to  a single  quadrant  of  the  complex  impedance  plane) 
is  proved  and  discussed. 

Using  methods  analogous  to  those  developed  for  the  corresponding  waveguide  problems 
by  J.  Schwinger,  the  lumped  circuit  representing  the  change  of  cross  section  ( and  of 
material)  in  circular  rods  is  shown  to  be  a simple  series  impedance  and  a variational 
expression  is  obtained  for  it  in  tonus  of  ths>  higher  mode  amplitudes  of  the  temperature 
field  in  the  circular  junction  region.  A set  of  linear  equations  gives  these  higher 
mode  amplitudes  where  the  coefficients  are  Bessel  function  suras.  Numerical  procedures 
for  evaluating  these  Bessel  function  sums,  and  in  particular  the  asymptotic  forms  of 
these  sums,  are  considered  in  detail.  Numerical  values  of  the  impedance  are  explicitly 
calculated  for  a range  of  cross  sections  and  frequencies  but  no  change  of  material; 
but  data  is  tabulated  from  which  the  impedance  for  a change  of  material  as  well  as  a 
change  of  cross  section  may  be  obtained  by  R simple  additional  calculation. 

The  static  heat  flow  problem  down  a rod  with  change  of  cross  section  and  material 
is  shown  to  be  solved  the  same  limped  series  impedance  as  the  dynamic  problem,  but 
evaluated  at  zero  frequency. 

Analysis  is  made  of  the  temperature  field  of  tho  ring  source  which  corresponds 
to  the  heat  produced  in  a thin  resistor  wrapped  around  the  surface  of  a rod. 

A tabulation  of  symbols  and  of  mathematical  information  relating  to  the  numerical 
calculation  is  given  in  the  appendices. 
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INTRODUCTION 


Recent  work  on  a method  of  measuring  thermal  properties  of  metals  which  uses 
non-statiouary  temperature  fields,  or,  mere  precisely,  which  usee  temperature  waves 
in  cylindrical  rods,  has  given  rise  to  questions  of  the  reflection  of  such  wavo^  at 
various  terminations  or  discontinuities . (1)  This  report  discusses  the  calculation  of 
such  reflections,  as  well  as  the  general  physical  background  of  temperature  waves, 
and  carries  through  a calculation  of  the  reflection  at  a change  of  cross  section 
between  two  circular  cylindrical  rods. 

In  Section  2 the  equations  of  the  temperature  field  TC^.tO  are  set  up  on  the 
basis  of  explicit  physical  assumptions,  and  the  characteristic  solutions  or  modes 
obtained  in  circular  cylindrical  uniform  rods  with  free  sides,  for  harmonic  tine 
dependence.  A table  is  given  of  values  of  diffusiv'.ties  and  other  thermal  properties 
at  various  temperatures. 

Si  Section  3 an  analysis  is  made  of  the  behavior  of  TU,t)  with  various 
boundary  conditions  and  heat  sources,  both  stationary  and  time  varying;  only  the 
form  of  T(f,t)  after  the  transients  are  over  is  considered,  and  for  simplicity 
only  the  principal  mode  is  used.  The  (harmonically)  time-varying  problems  are 
expressed  as  a superposition  of  a static  problem  satisfying  Laplace's  equation, 
and  a dynamic  problem  satisfying  the  wave  equation,  but  with  simpler  boundary 
conditions . 

In  Section  U the  transmission  line  analogy  is  developed  for  the  description 
of  the  behavior  of  the  principal  mode,  and  equivalent  voltages  and  currents,  chosen 
as  analogous  respectively  to  the  temperature  and  the  heat  flew,  are  introduced. 

The  effects  of  discontinuities  are  rigorously  represented  as  lumped  circuits  in- 
serted in  the  transmission  lines.  A condition  chi  the  characteristic  impedances  is 
obtained,  which  must  lie  satisfied  if  reciprocity  is  to  be  shown  by  the  linear  rela- 
tions between  currents  and  voltages  at  the  discontinuity,  hence  is  required  for  a 
circuit  representation  of  those  linear  relations.  The  transmission  line  formulas 
are  given  for  the  calculation  of  the  reflection  coefficients  due  to  various  inserted 
lumped  circuits  and  terminations  of  the  line.  A general  discussion  and  proof  are 
given  of  tne  theorem  ttmt  all  thermal  impedances  lie  in  the  first  quadrant  of  the 
complex  impedance  plane,  hence  are  resistive  and  capacitative.  Thermal  impedance 
here  is  the  ratio  of  the  (dynamic,  complex)  temperature  to  the  (dynamic,  complex) 
heat  flcr  through  a cross  section, of  a principal  mode  thermal  wave,  and  corresponds 
to  the  voltage  to  current  ratio  of  the  analogous  transmission  line. 

In  Section  5 the  mathematical  analysis  of  the  equivalent  circuit  for  a change 
of  both  cross  section  and  of  material  at  a junction  of  two  circular  cylinders,  is 
set  up  and  carried  through.  The  characteristic  impedances  are  uniquely  fixed  by  the 
requirements  of  reciprocity  at:d  simplicity  of  the  equivalent  circuit.  An  explicit 
variational  expression  for  the  circuit  element  (which  is  just  a series  element 
between  the  transmission  lines)  is  obtained  and  expressed  as  a sum  of  higher  mode 
amplitudes.  Linear  equations  for  the  higher  mode  amplitudes  are  given,  with  co- 
efficients, which  are  explicit  Bassel  function  sums. 


(1) 


E.  Mendoza.  - Bull.  Am.  Fhys.  Soc.  28  18  (Jan.  22,  1953) 

II.  Waldron  and  M. A. Merlin  - Bull.  Am.  Phys.  Soc.  28  23  (April  30,  1953) 


In  Section  6 , the  various  formulas  are  collected  and  put  in  a form  suitable 
for  numerical  calculation.  In  particular,  the  asymptotic  forms  of  the  Bessel  func- 
tion sums  for  the  CL<n  are  evaluated  and  used  to  simplify  the  calculation.  Calculated 
values  of  the  circuit  element  Z and  reflection  coefficient  are  given  for  various 
selected  values  of  the  physical  parameters.  Numerical  data  are  given  permitting  u imp- 
calculation  of  Z for  all  changes  of  cross  section  and  of  material  and  for  all  fre- 
quencies of  practical  interest*  The  general  coefficient  C^m  for  not  equal  to  m. 
is  shown  to  depend  simply  on  and  CorA,  so  that  the  labor  of  computing  the  square 

array  of  coefficients  is  much  reduced.  The  functions  used  by  Hahn  are  explicitly 
evaluated  in  terms  of  the  earns  functions  used  for  the  discussion  of  the  asymptotic 
forms  of  Bessel  function  sums. 

Section  ? discusses  the  static  problem  for  the  change  of  cross  section,  namely 
the  heat  flow  through  the  junction  vhen  the  ends  are  maintained  at  fined  (different) 
temperatures.  The  solution  is  shown  to  follow  immediately  from  the  same  equivalent 
circuit  representing  the  junction  as  in  the  dynamic  problem  (but  with  circuit  para- 
meters evaluated  at  zero  frequency). 

Section  3 derives  the  temperature  field  for  the  physically  realistic  ring  source 
(corresponding  to  a thin  electrical  resistor  around  the  rod)  which  produces  heat  uni- 
formly over  the  area  of  a ring  on  the  surface  of  the  rod.  The  solution  is  readily 
found  by  using  the  appropriate  Green's  function  in  the  cylindrical  region,  and  is 
expressed  as  mode  expansion.  Of  particular  interest-  are  the  phase  and  amplitude  of 
the  principal  mode  wave,  to  be  respectively  compared  with  the  phase  of  the  source  and 
the  amplitude  of  the  next  higher  mode. 

Appendix  I contains  a tabulation  of  the  symbols  used  in  this  report  and  their 
definitions,  listed  in  order  of  their*  first  appearance  in  the  body  of  the  report. 

Appendix  II  gives  procedures  for  evaluating  a definite  integral  which  appears  in 
the  sum  formulas  for  the  asymptotic  Bessel  function  seriee. 

Appendix  III  is  a sample  numerical  calculation  of  the  coefficient  C0I  using  the 
formulas  for  asymptotic  summation. 

Appendix  17  tabulates  various  numerical  quantities  useful  in  the  calculation  of 
the  lumped  series  impedance  Z or  the  reflection  coefficient  <R  for  the  change  of  cross- 
section.  The  numerical  values  of  Z and  (R  for  a range  of  values  of  the  ratio  of  radii 
and  a range  of  frequencies  (but  no  change  of  material)  are  tabulated. 
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EQUATIONS  AND  MODES  OF  TUB  TEMPERATURE  FMJ 

-A  / OufGS  - 

The  phenccnnolcgica^^guati.on  for  ’mat  flux  Q.  ( /car  sec)  introduces  the 


thermal  conductivity  K ( "7cm.sec.deg,) 

Q.--K7T 

while  conservation  of  (heat)  energy  requires 

-V-GL- 


g.i) 


where  -dai.sity  (E^/cm^),  C -specific  heat  ( /gm.°),  T-temperature  (deg.  Cent . ) , 

t«tin».  (secs. ) Putting  (2.1;  in  (2.2)  gives  the  basic  differential  equation  for  the 
temperature  field 


(Jl-3) 


which  will  always  be  treated  in  the  approximation  that  K is  independent  of  position 
(except  for  abrupt  discontinuities  at  junctions  of  materials  which  will  be  considered 
explicitly) > and  thus  Independent  of  ? as  well.  Then  (2.3)  takes  the  form 


/D){hT/$l)  ■ K/9C^  - /sec.; 


(?•*) 


The  solutions  of  (2.4)  for  T(X  X)  determined  by  appropriate  boundary  and  initial 
conditions,  are  the  quantities  of  interest.  In  particular  the  solutions  in  homogeneous 
cylindrical  rods  are  the  basic  mathematical  expressions  that  are  needed  in  the  problems 
to  be  considered.  Hence  we  give  new  the  "modes”  of  the  temperature  field  in  cylindrical 
rods  for  the  case  of  harmonic  time  dependence.  For  simplicity,  and  because  the  immedi- 
ate application  is  to  this  case,  we  discuss  circular  cross  sections,  although  a general 
theory  for  arbitrary  cross  sections  can  readily  be  given.  Also  for  these  reasons  we 
restrict  attention  to  the  circularly  symmetrical  modes  (independent  of  azimuth). 
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•V 


As  usual  we  introduce  the  complex  harmonically  varying  temperature 

Tlr.t)  - T“’(r)&iu,X 


(I- 51 


(where  ^ is  a general  position  vector  and  the  true  physical  temperature  is  either  the 
real  or  the  imaginary  part  of  T (r,t>  ),  and  introducing  cylindrical  coordinates 
v*  ind  % (2.4)  becomes 


} IT* 


M1 


v ia-r1. 

7 17  p ' 
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The  possible  solutions  of  (2.6)  ore  restricted  by  the  boundary  condition  on  ths 
cylindrical  surface,  at  radius  r«P..  Neglecting  radiation,  a very  good  approximation 


1 ! 


ft»n-«,TWrp  »«•»- 


for  the  lew  temperatures  (a  few  degrees  Kelvin)  in  mind,  this  condition  is  that 
Q/n  ■■heat  flow  normal  to  the  surfaced,  hence  from  C1*0 


vr 

a w 


(,a-7) 


I 
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(2.7)  will  then  restrict  the  possible  radial  mode  functions. 


Separation  of  variables  by  introducing  a product  function  for 
leads  immediately  to 


Tl\r)*  flr)  tU) 


^r) 


c. 


ico 

D 


U.9) 


U*  I'3) 


i 


The  solutions  of  (2.10)  are  the  Bessel  functions  of  zero  order,  ) » M<yC^T«  ^), 

The  latter  are  excluded  by  the  requirement  of  finiteness  at  r*0,  and  they* 
are  then  fixed  by  the  boundary  condition  at  r-R  giving 

i >0*0  , i.mft , 


(The  are  tabulated  in  Appendix  IV) 

Thus  the  temperature  modes  are  . v.  - . *. 

Tlr,t)»  * 


where 


t?.n) 


The  orthogonality  and  completeness  of  these  modes  will  be  introduced  laier, 
but  Just  now  we  note  that  in  the  z direction  they  take  the  form  of  propagating 


5 

(and  attenuating)  waves 0 In  particular,  the  lowest  mode  is 


T(y,t)*e  .e’jj'g  ’nr 


sO^iK 


U-B) 


Note  that  even  the  lowest  mode  attenuates,  although  less  strongly  than  higher  modes, 
and  that  it  is  uniform  over  the  cross  section,  thu3  obviously  satisfying  the  boundary 
condition o In  the  following  it  will  be  called  the  principal  mode  and  it  will  in  fact 
have  very  much  less  attenuation  than  the  next  higher  mode  and  all  others o Thus  its 
attenuation  is  a pure  constant  of  the  material,  whereas  for  all  higher  modes  the 
attenuation  is  increased  by  a geometrical  constant,  Yn  , arising  essentially  from 
the  difficulty  in  fitting  a long  (radial)  wavelength  into  a narrow  tube,  and  exactly 
analogous  to  the  "cutoff”  condition  for  electromagnetic  waves « It  is  the  vave  form  of 
the  modes  in  the  z coordinate  that  gives  rise  to  the  transmission  line  analogy  deve?.oped 
in  the  next  section. 

Finally  we  note  the  static  forms  of  the  temperature  modes,  obtained  by  letting  u>  , 
hence  also  'tf.  , approach  0.  (2.6)  becomes  Laplace’s  equation  and  (2.12)  (2.13)  take  the 

form 


TOr')  * 

Table  of  Thermal  Properties 


J.(y*r)£  r,  r\x  - 


(principal  static  mode) 

(higher  static  modes) 


U-M) 


On  the  following  pages  are  tabulated  values  of  the  diffusivitv,  thermal  conductivity, 
and  specific  heat  (per  unit  volume)  of  a number  of  metals  and  some  other  materials  at 
liquid  helium,  liquid  nitrogen  and  room  temperature,.  These  will  give  some  idea  of  the 
magnitudes  of  the  physical  quantities  in  the  above  discussion  and,  particularly  of  the 
diffusivity  '"or  various  materials  as  a function  of  temperature. 


The  wavelengths  A of  temperature  waves  follow  directly  from  the  values  of  diffus- 
ivity and  frequency.  Since  these  wavelengths  are  the  basic  measured  quantities  and 
essentially  determine  the  behavior  of  the  temperature  field,  they  are  also  tabulated  at 
a frequency  of  100  cycles/sec  with  a conversion  table  for  other  frequencies.  Since  the 
amplitude  of  a temperature  wave  drops  by  a factor  of  Q1**  or  about  0.2&  per  wave- 

length the  useful  length  for  measurement  is  les3  than  (A/a.)vhile  isolation  of  two  points 
from  each  other  is  effectively  achieved  at  about  a wavelength.  Hence  A must  be  neither 
too  long  nor  too  short  for  convenient  measurement. 


At  room  temperature  and  liquid  nitrogen  temperature  A is  a fraction  of  a centi- 
metre, which  is  uncomfortably  small  for  measurement.  However,  due  to  the  great  in- 
crease of  D at  helium  temperatures  (associated  with  the  decrease  in  heat  capacity  below 
the  Debye  temperature),  A,0,  has  the  convenient  size  of  a few  centimetres.  In  a few 
cases  of  very  pure  materials  A,,0  is  somewhat  Loo  large,  (e.g.  Sn,  Cu,  Mg,  Zn)  but 
increasing  the  freauency  can  easily  bring  these  down  to  a suitable  range.  Thus  wave- 
length considerations  restrict  the  temperature  wave  method  to  being  primarily  a low 
temperature  technique  at  present,  not  to  speak  of  other  technical  difficulties  connected 
with  temperature  measurement  and  behavior  of  transients  at  higher  temperatures. 


Values  of  Dif fusivities  and  Other  Thermal  Properties  at  Helium.  Nitrogen  (approx 


Room  Temperatures. * 

/C  « thermal  conductivity,  given  in  vstts/cm  deg 

C5  « heat  capacity  per  -unit  volume,  given  in  joule a/cs?  deg 

- /?  Cp  , /=>  ■ density  (gm/cm^)  5 * specific  heat  (jovQ.es/gm  deg) 

» diffusivity  , given  in  (cm  /sec) 

Au  - wavelength  of  temperature  wave  at  frequency  i>  (given  at  v - 100  cycles/sec) 


00  1000  i 2000 


; watts 
T (°K)  K\cm  de 


100 (cm 


nancy  or  sample 
for  measurements 


raphite 

iamond 


3.18  x 10"^ 
5.03  x 10-4 
7.54  x 10"4 
1.5$ 

2.42 


1.50 

2.00 


7.34  x 10~2 
1.77 


0.236  x lO4 
0.249  x 104 
0.212  x 104 
1.39 
0.827 

0.134 

9.72  x 10“2 
6.64  x lO"2 

0.665 

Or  464 

0.637 

3 .68  x 102 


7 

7 

6 

0 


0.32 

0.13 

0.11 

0.091 

0.29 

0.24 

0.28 


99.994  % 


0.41 

0.46 


0.152  x lO'2 


0.316  x 10 


watts  v 


K\cm  de 


II 


A 100  (cm) 


5.46  x 10“4  2,01  x 104  50 

9,40  x 10'^  1,49  x 104  43 

2,52  1,82  0,48 

3,42  - 1,11  0,37 


(normal) 


(super eon) 
ducting 


10,63  x 10"4  0 .,136  x 104 

0,292  x 10“J  7,53  x 102 

0,665  x 10“2  4,14  x 1C2 

10,63  x 10“4  7,05  x 102 

0,346  x 10*2  5,93  x 102 


(nonaaj) 


( super  con-) 
ducting 


(supereon-) 

ducting 


0,35 

0,5 

0,671 

13 

19 

17 

0.384 

0,346 


10,63 

x 

0,292 

x 10-- 

0,665 

x 10“2 

10,63 

x 10“4 

0,346 

x 10-2 

2,65 

0,141 

x 10“2 

0.218 

x 10“2 

0,300 

x 10”2 

3,51 

11,6S 

x 10~4 

a, 2 

x 10”4 

84,7 

1.32 

1,45 

x 10"4 

0,169 

x 10’2 

0,481 

x 10  “2 

1,10 

x 10”2 

2,215 

x 10“4 

3,89 

x 10“4 

5,99 

x 10”4 

0,222 


1,61  x 102 
1,67  x 102 
0,191 


1,12  x 104 
0.461  x 104 
0.201  x 104 
0,292 
0,238 

5,92  x 102 
7,27  x 102 
2,73  x 102 


37 

24 

16 

0,19 


99,99 

W 


99,95  % 


T (°K) 

/vatts  ' 
Kv,cm  deg. 

/Cm2\ 

D\»ec) 

A/(,j(cm) 

r.TTSJ 


t y of  aanple 
for  measurement 


(normal) 


( superoon-) 

ductir.g 


0.228 

x 10  * 

0.349 

x 10"2 

0.475 

*i  n-2 
X 

2.42 

3.95 

0.314 

x 10"2 

0,477 

x 10“2 

x 10* 
x 102 


99.997  % 

ft 


99.995  % 

n 


X 102 


18 

18 

17 

0,17 


99.999  % 

n 

n 


22 

17 

14 

0.51 

0.47 


99»99  % 


( supercon-) 
ducting 


(normal) 


99.98  % 


T (°K) 


9 


K\cm  de 


cisas*)!  D(Ss) 


nty  oi  sample 
for  measurement 
sf  X 


2 

36 
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3 

51 

white 

4 

55 
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0 ,799 

291 
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2 

22 
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3 

46 
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4 

55 

2 

0.015 

3 
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4 
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2 

0,01 

3 

0,015 

4 

0.02 
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1 .46 

2 

5 

3 
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4 

10 

118 

1,16 
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1.10 

2 
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3 

4.0 

4 

6.0 

2 

0.02 

3 

O.03 

4 

0.04 

1 

90 

0.434 

c 

lpt)2.l86 

17=2 

20^9 

4 
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3o24 

x 10  * 

7=40 

x iO"4 

15  = 5 

x 10"* 

1.23 

1.48 

3=24 

x 10"* 

15  c 5 x 10“4 

0.0669x  10“2 
0.118  x 10"2 
0.181  x 10“^ 

0a0669x  10~2 
0.118  x 10“2 
0.181  x 10~2 

2.74 

1.96  x IQ'1* 


e*  #«•/ 

t 0 AO 

x 10-4 

2.21 

2,77 

1.96 

x 10"^ 

4.15 

x 10~4 

7.76 

x 10-4 

0 

c 

0 

r- 

-0 

x KT2 

0,0822x  10~* 

11.1  x 10* 
6.39  x ICA 
3 = 56  x 104 
0.649 
0.435 

6.79  x 10^ 
10 

3=56  x 1C*1 


2.54  x 10^ 
1.81  x IQ1* 
1.29  x 10** 
0.524 
0.397 

1.02  x 10^ 


0.772  x 10 


118 

93 

67 

0.28 

0.23 


56 

48 

40 

0.26 

0.22 


99=996  % 

ft 


99=99  % 

ft 


99=9  % 

n 


99.997!^  Single 
*»  crystal 


~ 98.0  % 

ft 


0.249 

1.02 


Assumes  same  K 
as  for  He I 


* Data  on  heat  capacities  from  C.  A.  Shiffraan  ’’The  Heat  Capacities  of  the  Elements  Below 
Room  Temperature”,  G„  E.  Report,  October,  1952. 

Data  on  helium  temperature  thermal  conductivities  from  references  given  in  J.  Olsen 
and  H.  Rossn’oerg  '^Thermal  Conductivity  of  Metals  at  Low  Temperatures”,  Advances  in 
Physics  2,  28  (1953)= 

Other  data  from  the  Handbook  of  Chemistry  and  Physics. 


at  w^juga  raji&aa 


STATIC  AMD  DYNAMIC  PASTS 


Distributions  for  Static  Boundary  Conditions  and  Heat 


Sources. 


The  overall  purpose  of  Section  3 is  to  show  how  the  general  heat  conduction 
problem  with  time  varying  boundary  conditions,  separates  mathematically  into  a 
static  and  a dynamic  problem  for  each  of  tho  various  boundary  conditions.  Although 
the  measurements  are  primarily  concerned  with  the  dynamic  problem,  this  analysis  is 
necessary  because  the  two  problems  are  usually  mixed  in  the  physical  situation. 

Also  the  nature  of  the  solutions  is  not  always  obvious,  since  ixi  some  cases  there  is 
no  true  static  solution,  independent  everywhere  of  time.  The  dynamic  problem  never 
has  such  difficulties,  however. 

Accordingly  we  start  in  Subsection  (3a)  with  * tabulation  of  solutions  of  sous 
simple  static  problems  with  various  boundary  ccnditious.  These  will  be  useful  in 
discussing  the  more  complex  time  varying  problems,  and  will  develop  some  feeling  for 
the  physical  behavior  in  the  situations  of  interest. 

Two  types  of  physical  boundary  conditions  commonly  occur,  and  we  illustrate 
their  effects  ok.  the  temperature  field  by  typical  simple  (cna  dimensional)  situations. 
Either  T may  be  fixed  over  a given  surface,  say  by  means  of  a reservoir  (of  infinite 
heat  capacity  and  heat  conductivity) , or  CdT/dn)  may  be  fixed  over  the  surface. 

(A  still  more  general  boundary  condition  fixes  the  ratio  of  T to  (dT/dn).  This  is, 
in  fact,  the  radiation  condition  unto  a vacuum,  and,  as  we  shall  see  below,  amounts 
to  specification  of  the  "thermal  impedance"  of  the  boundary.  Some  consideration  is 
given  to  this  boundary  condition  later,  but  the  possible  physical  situations  are  adequately 
illustrated  by  +-he  special  cases  above.)  Two  possibilities  for  (dT/dn  ) are  conveni- 
ently distinguished  j either  (cfT/dn  ) ■ 0 , corresponding  to  an  isolated  or  free  surface 
across  which  no  heat  flows,  since  by  (2.1)  G I'Ti  *»0j  or  (ST/dn)  is  finite  and  heat 
does  cross  the  surface.  The  latter  situation  corresponds  to  a source,  it  >rill  be  treated 
by  explicitly  introducing  the  source,  and  in  fact  generalized  to  the  case  of  an  inter- 
nal source.  (A  third  type  of  termination  which  may  nominally  be  considered  a boundary 
condition,  is  infinite  extension  of  the  rod  in  so me  direction.  True  static  solutions 
for  this  case  do  not  exist,  but  for  fixed  after  long  enough  time,  static  bdiavior 
does  set  in.) 

Three  combinations  of  the  boundary  conditions  T fixed  or  ( oT/otO  • 0 are 
illustrated  in  the  figures  for  a section  of  uniform  rod.  These  are  all  static,  i.e. 
time  independent,  solutions  satisfying  crT/a-^  ■ 0 . 


fi<l  I.  * to/*.  \ 


J.  1 -C.tt.Vfc  A 


i,  3 - co.!,fc  3 


constant  value  depend*  on 
L ln\t 


V V Y 

In  Fig.  1 the  ends  are  at  fixed  values T,  and  77  , T.?Ta  , and  the  steady  state 

to  T-T-  cy/AKT.-Tj  ' 1 

In  Fig.  2 one  ouli  ~s  at-  "y  hue  other  end  is  free,  and  the  final  steady  state 
is  T » T4 . 

In  Fig.  3 both  ends  are  *Ve®,  and  any  constant  value  of  T is  a solution. 

(In  the  diagrams  the  representation  adopter  will  usa  solid  lines  for  free  walls,  and 
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dotted  lines  for  constant  temperature  surfaces,  a choice  that  seeros  natural  la  the 
sens©  that  free  walls  are  adiabatic,  whereas  heat  passes  across  the  constant  T 
surfaces.  The  logical  extension  of  this  notation  will  use  a double  solid  line  for 
a heat  source,  which  corresponds  to  a fixed  value  of  (dT/in  ),  or  of 
the  discontinuity  in  ( for  an  interne!  source,  a wavy  solid  line  far  a time 
varying  source  efci.  These  extensions  will  be  illustrated  later.) 

To  complete  the  tabulation  of  static  problems,  consider  the  steady  state 
temperature  distributions  whan  a source,  Q,  , is  introduced  at  into  the  rod 
(of  cross  section  ares  4r.  For  this  purpose  to  use  a simplified  source  which 
introduces  heat  uniformly  over  the  cross  section  (i.e.  a circular  disc  source  of 
zero  thickness.)  With  such  a source,  T remains  uniform  over  the  cross  section,  but 
has  a discontinuous  derivative  at  given  by: 


C3-0 


The  three  previous  cases  become  those  in  Figs.  4,5,6.  Thus  in  Fig.  4 T has  two 
linear  sections  which  intersect  at  Q,  and  are  explicit!  giveyiby  (3»2) . 
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Finally,  in  the  situation  in  Fig.  6 no  true  steady  state  exists  since  there  is  no 
way  to  dispose  of  the  heat  without  heating  up  the  rod.  However,  the  transients 
depending  on  initial  conditions  do  disappear  and  the  solution  approaches  a simple 
form  made  up  at  any  instant  of  two  parabolic  sections  intersecting  at  Q.( 
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As  time  passes  all  points  of  the  rod  Increase  their  T at  ths  same  rate,  so  that 
this  situation  might  be  characterized  as  a "quasi  tatic"  state.  It  is  not  time 
independent  as  in  previous  cases,  but  at  least  is  j hanging  uniformly  (not  periodically) 
with  time,  and  is  independent  of  ths  initial  conditions  (except  for  the  constant  C ). 

Note  that  (3.3)  and  (3.4)  provide  the  solutions  to  the  two  problems  with  a source 
at  ths  boundary,  since  ths  source  may  move  to  the  free  end  without  affecting  the 
solution,  e.g. 


(.3.5) 


T-  v ftm 


The  various  problems  with  semi-infinite  termination  of  the  rods  are  solved  later. 

Periodic  Temperature  Distribution  far  Harmonically  Varying  Sources  and  Boundary 
Conditions.  ~~  ” ” 

Now  suppose  the  source  has  a harmonic  component  so  that  the  heat  input  is 
k e *L°V.  The  solution  is  obtained  as  the  superposition  of  solutions  of 
v ° two  problems  which  may  be  designated  as  the  static  and  the  dynamic 

problems.  The  static  problem  has  the  steady  state  solution  set  up  by  &,  and 
already  given.  The  dynamic  problem  requires  solution  of  a pure  wave  propagation 
problem  set  up  by  a source  Q^e*^6  but  with  the  simpler  boundary  condition  T»0 
on  any  surface  held  at  fixed  T.  This  has  a solution  with  harmonic  time  variation 
(after  the  -transient  has  died  out),  which  is  independent  of  the  initial  conditions, 
and  is  conveniently  call  the  periodic  state  (a  suggestion  due  to  B.  Mendoza). 

We  may  represent  these  problems  dia grammatically,  as  follows,  with  the  boundary 
condition  of  the  three  cases  corresponding  to  the  constant  source  problems  above. 
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The  general  solution  of  Case  (1)  now  has  the  form 

x»  rJ)  * tij'  aioJt 


v^ere  'P' and  T^are  solutions  of  the  static  and  dynemic  problems  (3.6)  and  hairs  the 
forms 
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Note  that  T^ln  (3.7)  is  Identical  with  (3.2). 

Similarly  Case  2 separates  into  the  static  problem  with  solution  in  (3.3) 
and  a dynamic  problem  with  the  solution: 
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Finally,  Case  (3)  separates  into  the  quasi-static  problem  (3.4)  and  the* 
dynamic  problem: 
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The  problems  with  time  varying  sources  solved  above  describe  the  essential 
physical  situation  in  the  heat  conduction  resulting  from  resistor  sources  fed  by 
A.C.  Although  a simplified  source  has  been  used  wriich  gives  rise  only  to 
principal  mode  waves,  less  symmetrical  sources  will  behave  similarly.  They  will 
yield  in  general  both  principal  and  higher  mode  waves  but  the  latter  may  be  ignored 
at  sufficient  distance  from  the  source  in  comparison  with  the  principal  mode  because 
of  their  heavier  attenuation.  {Sea  Section  8 for  an  analysis  of  the  ring  source. ) 

Another  and  simpler  group  of  nan-stationary  problems  which  are  not  quite  as 
easily  realized  physical^-  but  still  may  be  approached  in  some  circumstances,  assumes 
the  temperature  of  the  boundary  has  a harmonic  part  as  veil  as  e constant  part. 
Again,  these  problems  may  be  reived  siupiy  as  the  superposition  of  a static  and  a 
dynamic  problem.  Only  Cases  (1)  and  (2)  are  relevant  now  (Case  (3)  has  no  surface 
at  constant  T ) Soh.tionm 
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ri”"Uy,  it  *111  Lc  useful  to  note  trie  solutions  to  the  above  source  prcllaaa 
w3.t.;  iu: other  termination  - namely  with  infinite  extension  of  the  conductor  in  a 
dire  away  .from  the  source  (indicated  by  an  arrow  and  open  end).  There  are  four 

such  problems: 
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In  the  static  part  of  all  these  problems  with  semi-infinite  terminations,  the 
transient  dies  away  in  a time  depending  on  the  distance  from  the  source  or  temperature 
reservoir*.  Thus  the  solutions  depend  on  the  error  function  of  the  argument  (y-/^ 
where  c-.'f  x - < <±/nj \ ">A/  , which  rapidly  .approaches  its  limiting 

value  of  / when  the  argument  exceeds  ! . Accordingly 

if  at  fixed  y the  temperature  approaches  a constant  (time  independent)  value  for  long 
enough  time  (the  time  will  be  proportional  to  ),  we  have  called  the  solution  static 
(perhaps  asymptotically  static  would  be  better).  If,  however,  7~ "''approaches  a value 
rising  steadily  with  the  same  rate  at  each  point,  as  ±r.  (3*16)  arid  (3.18),  the  solu- 
tion has  been  called  quasi-static,  in  analogy  with  (3.4)= 


These  examples  illustrate  that  various  heat  flow  problems  with  time  varying 
sources  can  be  solved  as  the  superposition  of  two  problems.  One  is  a dynamic  problem 
with  simplified  boundary  conditions  - cither  7" or  in  which  the  periodic 

solutions  are  made  up  of  temperature  waves  drive  . by  the  periodic  portion  of 

the  source.  The  other  is  a static  problem  which  keeps  the  more  general  boundary 
condition  with  respect  to  T , In  certain  cases  no  true  static  solution  'sts  far 
this  problem,  but  then  a quasi-  stati.c  state  in  which  the  T of  all  points  a s at 
the  same  rate  does  exist.  Ho  such  difficulties  arise  in  the  dynamic  problem  In 
all  these  cases,  the  transient  solutions  which  depend  on  tteaxact  initial  con_  .tions, 
are  not  of  interest  here „ 


(3c)  Classification  of  Problems  Treated  and  Their  Equi  valent  Simpler  Problems, 


Although  all  the  discussion  lias  been  given  with  simplified  sources,  uniform 
over  the  cross  section,  this  restriction  does  not  affect  the  essential  nature  of  the 
solutions  under  the  various  conditions  considered.  The  examples  discussed  above 
should  be  typical  cf  all  physical  situations  of  heat  conduction  in  cylinders  with 
the  boundary  surfaces  at  fixed  7~  , at  fixed  flux  (i.e.  fixed  value  of  ^T/in  ),  fbee 
(dTAn  -&)  t or  with  time  varying  7",  or  time  varying  flux.  In  addition  the  boundary 
condition  of  semi-infinite  extension,  and  the  more  general  case  of  internal  sources 
have  been  considered. 


The  problems  treated  may  be  systematically  classified  as  follows,  where  the 
solution  to  the  complete  time  varying  problem  on  the  left  contains  all  the  simpler 
problems  on  the  right.  Brackets  join  the  two  problems  (static  and  dynamic)  whose 
superposition  solves  the  time  varying  problem.  The  static  problems  include  internal 
source  problems  with  all  six  combinations  of  the  three  boundary  conditions,  fixed  77 
free,  or  semi- inf inite  (/(.,’ . . 17  *-*X  In  case  the  internal  source 

occurs  with  a free  boundary,  these  may  be  com- 

bined to  give  an  end  source  with  the  three  boundary  conditions.  When  the  sources  go 
to  zero,  the  six  possible  distributions  with  no  sources  are  obtained.  Finally,  the 
time  varying  boundary  TJs  give  on  the  one  hand  three  static  problems  which  are 

included  in  the  six  no  source  ca3e3  just  mentioned,  and  have  fixed  7~  at  one  end  and 
three  boundary  conditions  at  the  other,  and, on  the  other  hand,  three  simple  dynamic 
problems,  classified  under  end  source  proble.no,  since  the  disturbance  is  produced 
at  the  variable  7~  end.  The  notation  used  in  the  diagrams  is  the  same  as  that 
introduced  above  and  ic  defined  again  in  Appendix  1, 
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THE  TRANSMISSION  LINE  ANALOGY.  EQUIVALENT  CIRCUITS  XP.  DISCONTINUITIES 
AMD  THE  CALCULATION  CF  REFLECTIONS. 


(4a)  Transmission  Lice  Equations  and  the  Principal  Mode  Temperature  Wave 

The  form  of  ths  z dependence  of  any  of  the  temperature  modes  as 
either  running  or  standing  waves  permits  immediate  representation  of 
these  modes  by  a transmission  line.  The  usefulness  of  this  procedure 
rests  on  the  practical  possibility  of  considering  only  the  principal 
mods  and  its  transmission  line,  since  the  greater  attenuation  of  the 
higher  modes  damps  them  out  relatively  quickly.  Thus  communication 
between  discontinuities  such  as  sources,  junctions,  or  detectors  can 
be  considered  to  take  place  entirely  via  the  principal  mode  provided 
these  are  spaced  sufficiently  far  apart.  The  effects  of  the  discon- 
tinuities can  in  fact  be  represented  by  a lumped  element  circuit  which 
simply  describes  by  means  of  linear  relations  the  net  effect  of  the 
discontinuity  on  the  lowest  mode  voltages  and  currents  appearing  at  the 
circuit.. 

We  consider  now  the  principal  mode  and  describe  it3  behavior  by 
means  of  a transmission  line  - other  modes  may  be  treated  similarly. 

In  general,  (omitting  the  factor  e'/u>,::  and  dropping  superscript  ) 


stT  * Jt  X(A  e 1 *>)  - XUtfcos  X?  -«  s>nK$. ) 

a* 

c\  = A r 3 , jS  -£ 


to  be  compared  with  transmission  line  equations  and  solutions 

* / r z,l(f)  ; Vfi )’  K'osXj.**. 4 /.  s/nX? 
tUJU  .iXK  v(2)  , If 2)  = I cos  if./  *4.  X V.  sir,*? 
V.*V(o)i  Z„  s l(o)  j 


(4-0 


(4.2) 


Two  j.  i.Cvj-oi.Vr>iV3  of  V($)( and  I( $)J  are  possible, namely 

with  T ($■)  or  with  ( r-nd,  correspondingly  [(#)  w/th  , or-  T ) 

There  seems  to  La  no  important  reason  to  prefer  one,  although  we  note 
that  (9  '/c?  has  some  analogy  to  the  transverse  electric  field  in 
waveguide,  ‘since  both  vanish  at  the  free  surface  of  a junction  with 
a cylinder  of  s^j.lor  cross  section  (on  the  metal  wall  in  the  waveguide 
o..io,  ar.a  or  the  free  surface  in  the  present  case).  Since  voltage 

juii  i.  ..otl-xC  i.Lfc.'.Cj 
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ax»j  related,  this  leads  to  identification  of  (a7/s^.)  with  \/(?)  ■ 

However,  a different  analogy  comjjares  , which  ia  proportional 

to  the  heat  flux,  with  the  electric  current.  Thus  a free  surface  Is 
now  an  open  circuit  and  no  heat  flux  ga  ts  through  it„  Because  there 
is  seme  precedentfor  the  latter  choice  (private  communication  from  E. 

Mendoza)  we  take  it  as  the  basis  of  the  analogy.  Another  advantage  is 
that  this  will  lead  to  the  natural  representation  of  a heat  cry:  city 
by  a capacitative  (and  resistive)  circuit  element. 

Thus  we  shall  put 

V(%  )=CT  * CY«  cos  jc  ? sm  X?)  (4'  j) 


where  proportionality  constant  C ia  introduced.  Comparing  (4«,3)  and 
(4„2)  yields  then 

J(j)  Z7X*:  ^ * -iCKUtfcosX.z-'xsinXz)  (4-4) 


where  two  constants  remain  to  be  fixed,  C and  YQ  . These  will  be  chosen 
later  to  simplify  ohe  equivalent  circuit  of  the ^dicontinuity  ( and  to  satisfy 
reciprocity,)  of  course  ia  already  fixed  and  oc  and  /S  give  the 

wave  amplitudes  , 

(4b)  Effects  of  Diaconuinuities 


To  illustrate  these  equations  and  to  lead  up  to  the  description  of 
a junction  between  different  rods, consider  the  effect  of  terminating  a 
simple  rod  in  various  ways  Thus  the  termination  may  be  T-0  (the  complete 
problem  including  the  static  part,  will  have  T-constant  )hence  V-0  and 
this  should  be  considered  a short  circuit.  This  gives  (talcing  the  origin 
at  the  junction) 

e •*,  - o 

Vfy)  - C *./9s 

! I(j)  - c Yl3c  oz  X £ 

y-° 


The  reflection  coefficient  (for  the  voltage  waves)  is  /p's  ~ yjociyi))*  ~f. 

Similarly,  if  the  termination  13  a free  end,  T/z  o. ) hence  L*0  and  the 

rod  may  be  considered  terminated  in  an  open  circuit  corresponding  to 

Then  the  voltage  reflection  is  R - = /*/  . Fin  illy,  if  the  line 

continues  without  interruption,  i.  e.  is  matched  at  any  pojLnt,  then  ii*G, 

hence  R-Oy  ^**3,  - C'*  , /(?)*  C ^ 

and  the  .line  ne.y  bo  considered  terminated  in  its  characteristic  impedance. 


In  thy  general,  case,  however,  in  which  two  rod 3 or  cylinders 
meet  in  a junction  of  arbitrary  shape  and  may  in  addition  have  differ- 
ent sizes  and  materials,  ( as  in  Fig.  1)  the  situation  may  be  represented 
by  a lumped  circuit  between  bwotransmission  line3  with  different  propa- 
gation constants  and  characteristic  impedances, (shown  in  Fig.  2) 


rod  / 


\\\ ' / 

•.  i \ .i  / 


. < . / / / / / i 

/,'//•  if i // // ■ 

•://  //////// 


■rod  R 


F/g  / Change  of  cross 
section  end  fabr  r/c  t) 
change  of  material  or  a 
juns nan  between  rods. 


7'" 0.1  J//')!  $.  . /<>// 


/in*-  / 


Is  A. 


in  ~ s / u equivalent 

. »j».  > / 

— tronsmissior  /ines  and  hn, pea 
w m r CitC'Jif  t or  the  junction 

between  rods. 


To  justify  the  use  of  the  lumped  circuit,  we  note  that  w©  can  write 
linear  relations  among  the  principal  mode  amplitudes  at  certain  reference 
planes  in  the  two  rods,  in  the  form 


v.-- 

V ‘ h.T.  rZnTJ 


[4 5) 


(in  14.5)  the  positive  directions  of  voltage  and  current  in  the  tve 
lines  are  taken  as  in  Fig  2.  i.e.  positive  current  flow  toward  increasing 
d in  both  lines,  positive  voltage  when  the  upper  conductor  has  the 
higher  potential.  Quantities  referring  to  rod  II  will  be  designated  by 
primes. ) Such  a linear  relationship  roust  hold  because  of  the  linearity  of 
the  entire  problem  and  the  fact  that  two  quantities  y 7,  IS  cr 
completely  determine  the  field  distribution.  To  construct  an  equivalent 
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lumped  circuit-  to  represent  tuis  relation  requires  in  addition  that 
a reciprocity  condition  be  satisfied,  ■which  here  takes  the  form 

- ~Z2,  fsze  4/4). 

(4c)  Derivation  of  Reciprocity  Condition 

The  satisfaction  of  the  reciprocity  condition  is  possible  in  general 
if  the  characteristic  impedance  satisfies  a particular  condition,  which 
may  be  found  by  a procedure  like  that  used  in  the  electromagnetic  case* 
Thus,  Green1 s theorem,  applied  to  any  volume  V of  the  conducting  rods 
with  surfaced  gives  for  arbitrary  functions  fa  on-J  , 


J(fav!<Pb-<Pb 


\/ifa)d\ '*/(<& 

JT 


(4.6) 


Put  for  <Po/  -p6  >0  ( 4-  ■-•)  two  different  solutions  to  the  dynamic 

problem,  Ta  and  Tb,  whic'n  then  satisfy  (% 7 i Xz)  7k&  *0 . 

Thus  the  vclu-re  Integra1,  in  (4.6)  vanishes,  and  since  ( J ) - O 

on  the  side  walls,  uhe  only  contribution  to  the  surface  integral  comes  from 
the  two  cross  sections*  Let  these  be  at  &ncf  5 = Z2  > L,  ; 

and  substitute  from  (4.  2).  (4  3),  (4  4), 


T--  ( ~J  COi r Xj.  t (i  Zo  T3/c  ) S>n  // 

(ij)  -- 


(4.7) 


(4=7)  gives  the  lowest  mode  forms  of  T and  which  are  dominant 

sufficiently  far  from  the  junction  , hence  may  be  assumed  to  hold  at 
3 ~ £-4  j Z 2 * Here  V0  and  I0  are  the  values  of  V find  I at  the 

reference  plane  z-0,  Then,  substituting  (4.7)  in  (4.6), 

/[(  c CS  XL,  7 i if.  Is?  s in  XL, )(-  Y~3Yz3  c /n  ML,  + (LKJLl  (s'co:  X L.) 

L < c c n C c 

..  ^ YsS cs  XL,  t ZYLIzP  sin/L ,J(  S/n  XL,  4 sLJpIc?  cos  X Lj/ 

- cirn/hr  expression  :jT  3.  *L2. 

n t~  ) ( 9 T 1 

where  ( / ~ ~~\  0 3)  L has  been  used,  and  4is  the  cross  section 

area  of  rod  I,  Now  the  left,  side  of  (4.8)  reduces  to 

A tprt1™  ^ *■*>*  v*aj 


A t~°  “ ' - 


4 ’ - -V  ” ■(■  vt>  ",~£o.h  l Oof 


(4.9) 
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where  the  linear  ecuations(4*5)  between  V and  I at  s&bO  in  the  two  rods 
have  been  used  . If  there  is  no  change  of  material  between  the  two  rods, 
but  just  a junction  region  and  possible  change  of  cross-section,  then 
Green’s  theorem  applies  at  any  two  cross-sections  Li  ,and  I»2  , while  (4*9) 
holds  sufficiently  far  from  the  junction  region.  Hence  (4*9)  would  give 
(note  that  73,  (*  + ) are  expressed  in  terms  of 

and  t2  ‘oc’  in  the  left  and  right  lines  respectively) 


A Hi  +Z*  Uo]--A~p‘L-ti  Ki+Jos  l/JJ  (4-./OJ 


or,  as  in  (4*9)  on  using  (4*5) 

4 }(  j f t t*  ' r -r  ' 7 /*  *7^  A 


Airrz'*[- z2/[-zJi„  +/„  ’i„]  (+.//J 


Hence 


jj  ^ a‘  z*  . , , . , 

YF77*  w">  9/ve  reciprocity 


relation  2*/2  * ~ i?2/  , which  is  necessary  if  a lumped  circuit  repre- 

sentation is  to  be  used. 

However,  if  the  material  changes,  this  proof  breaks  down,  since  the 
wave  equation  ( )t  *0  does  not  hold  in  the  regj^n 

where  the  material  changes.  The  proof  can  be  modified  if  the  change  it 
abrupt,  taking  place  at  a particular  cross-section  £ » ,> 

where  there  may  also  be  an  abrupt  change  of  section;  then  we  can  write 


r(2T)  -*'(£T) 

where  use  has  been  raae  ef^'Uie  continuity  of  heat  flow  * \<*  j- ^ 


as  in  (5*5)  and  that 
one  rod  ends  at  the  junction. 


■0  on  the  extra  (free)  walls  where 


Thus  the  general  condition  which  gives  reciprocity  even  for  a change  of 
material,  provided  this  happens  abruptly  at  a cross-section,  is 


f< A X7-  ..  /C'A'H‘Z.  ' 
C2  = iC')2 


(4.  /3 


Care  will  be  taken  to  satisfy  this  condition  later  when  .,he  change  of 
cross-section  in  circular  guide  is  explicitly  solved.  Note,  however, 
that  since  ) C is  complex  for  temperature  waves,  corresponding  to 

attenuation  of  the  principal  mode,  the  characteristic  impedance  will  also 
be  complex, 

( 4d ) The  Equivalent  Circuit  and  the  Calculation  of  Reflection  Coefficients 


Assuming  reciprocity  we  have 

V.  •-  Z.,7*  ' 

V.’  --2,. !4 


with  the  circuit  representation 


(4  14) 


L \ , L,‘ 

2.,  ~2U  rf,. 
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, -^VVytv^V— 

N l > 

Vo 
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^If  the  positive  direction  of  current  flow  in  the  primed  line  were  to  the 
left,  i.e.  if  we  int.rod;.  e the  current 

/v _ -cl  ±T  ^ c.‘  _ ±r 

iX'Zo  3 r j y where  ? ' r ~ y 

so  that  the  ‘ coordinate  has  the  same  relation  to  the  junction  as  s 
does,  then  the  sign  on  the  right  of  (&.  H)  or  (<4./2) 
would  be  reversed  (on  using  I0"  ) and  assuming  (4,13),  we  would  have 
reciprocity  expressed  by  Z,z  - JPi,  . 

Hence,  circuit  equations  and  equivalent  circuit  would  have  the  more 
symmetrical  forms 

= Z/)  I,  rZsi  'Fo  " 


(*  /r  j 


with  equivalent 
lo 

circuit 
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where  VZsVZ,  L"s  ~L\  Z.'.'-Z.,,  Z’h 
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~x  ■ 


£>?  yic-4-B)  xc .'  --*}) 


(4.  /9J 


where  (4.3),  (4.4)  and  (4.1)  have  been  u?ed=  (4.17),  (4.18) 

and  (4.19)  give  the  general  formula  for  the  reflection  coefficient 

J±K  - TTzlt  -f~  faz/z*) 

/-<%_  7*22  ~ £ IzL Zj  j-an^t  '7.)  . (4-.2&J 

(t~A  Yo'  Z't  ton*  I) 

Certain  special  cases  i (4.20)  are  of  interest , namely 

a)  Short  circuit  - V‘(4)~o  f/t*  7~(^ ) = O) , -O, 

or><J  (4.20)  jbecor-n? '.s 

/+<R  . y. 

/-*  ~ Z Z(Zz-+<Z •*,»*.'£)  (4.2!) 

b)  Open  ci  cuit  - I’tfJ'C  0 * (*Tf'&})  ~o).  ^ 


/+  fj 

/-  ft 


ir 


Z(Z,  -aZ'co  + K'Z) 


(4.22) 


c)  latched  termination  ~ (/€*  fOT /Tl'YiQtfKj  111  'ittLif’OCT&riStlC 

/mp^oco  'r  s$rni- inf injre  extension  of  /A-  /in?). 


?;*7‘ 

' ' — —t  —n 


It  ft  _ 2^  r Z,Z 

7-ft  2*  Z0{Zit~Z0) 

(4®)  General  Significance  of  the  Source  Problems 


(423) 


Finally,  wo  note  that  the  various  heat  conduction  problems  solved 
in  Section  3 which  use  simplified  sources,  and  uniform  cylinders  with 
the  three  possible  terminations  (constant  7"  = O,  semi-infinite ) 

may  be  given  a more  general  significance.  Both  discontinuities  in  ^cylinder 


and  general  sources  be  present  provided  eomiauiication  between  these 
and  the  boundary  crods  sections  takes  place  only  via  the  principal  mode 
An  example  provided  in  the  figure  is 


General 
tune  varying 
source 


Genera! 


discontinuity 


Static  Problem 


Dynamic  Problem 


Thus  the  general  problem  above  separates  into  a dynamic  problem  and 
a static  problem,  The  former  consists  of  the  end  conditions  on  the  two 
transmission  lines  (short  circuit  left,  open  circuit  right)  the  Iciest 
mode  part  of  the  source,  and  the  equivalent  circuit  of  the  discontinuity,  - 
The_  static  problem  consists  of  the  sra.tic  boundary  conditions 
( / - con  star/ U y * T'dj-  )*  the  steady  part  of  the 

source,  and  the  general  discontinuity.  Also,  strictly,  there  remains  the 
higher  mode  part  of  the  source  (time  varying)  which  simply  dies  away  without 
reaching  the  di:: continuity  - hence  without  generating  any  more  lowest 
mode  contribution.  The  physio  1 character  of  all  these  problems  is  included 
in  the  e:amples  in  Section  3*  although  it  may  not  be  possible  to  write  out 
the  solutions  explicitly. 

(4f)  The  Capacitative  and  Resistive  Character  of  Thermal  Impedances 


The  characteristic  impedance  and  the  constant  C in  (U°3  ) are 
fixed  later  in  (5d)  by  the  requirements  that  reciprocity  be  satisf  ed 
and  that  the  circuit  of  the  change  of  cross-section  simplifies  to  a 
single  series  element.  The  result  of  these  requirements  is  given  in 
(5.15)  and  makes 


(r  Ar  ) 


c =■•/ 


(■?  I fj 


with  the  natural  identifies  ion  of  l(J)  with  the  total  heat  current  down 
the  rod , Q(3-)j  and  V(3)  with  the  principal  mode  part  of  ~T(UO 

7S>(?)  ( c.f ■ n./6).  We  shall  anticipate  this  identification 

new,  in  order  to  set  down  in  definite  form  some  results  of  general  character 
on  the  behavior  of  the  impedance  seen  by  a temperature  wave  or,  briefly, 
of  a thermal  impedance 


From  (4»<24 ),  ptu  ting  in  the  complex  value  of  >£ 

<r.  ~ ( /+A)\f~-r  h.  (z+nk 

v c L> 


from  (2.13),  namely 

(4-mj 


where  (4»25)  defines  the  real  wave  number  A,  we  have 


2?  - 

-t-o  , 


KAA 


(4.26) 


(4.26)  shows  that  2T.  for  an y rod,  of  arbitrary  thermal  conductivity, 
cress  section  area,  or  wave  number,  i3  in  the  first  quadrant  of  the 
conplex  impedance  plane,  on  the  45*  lire,  i.  e.  thac  is  resistive  and 
capaciiacive . ( note  that  the  choice  of  the  time  factor 
makes  capacities  positive  imaginary,  the  opposite  of  the  usual  convention 
for  electrical  engineering,  but  more  convenient  far  problems  pf  wave 
propagation.)  Tills  result  is  readily  understandable  physically  since 
it  means  that  the  current  or  heat  flux  builds  up  before  the  voltage  <«• 
temperature  (remember  the  negative  time  exponent)  just  as  in  a condenser 
the  current  leads  the  voltage  * Both  the  heat  capacity  and  the  condenser 
are  reservoirs  to  be  filled  \rith  heat  energy  or  charge  respectively  before 
the  corresponding  "potential’*  appears. 


A more  general  result  follows  by  considering  the  class  of  impedances 
produced  by  transformation  of  some  ZJ  down  a length  of  thermal,  trans- 
mission line,  which  may  have  different  Z ona  Z-  From  (4.16),  this 
gives  rise  to  impedance  Zj 

n — - fonhiXX 

^ _ Zo  Zc  -Ton  <K.(c  _ £* 

/“‘A  K Do* A°n  X £ / - Os) # 
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In  (4.2?)  \~°  / 2*)  is  raal  arm  positive.  For  the  case  <?{?  ~ /, 

■fanh (2*/?„)  is  also  real  and  positive.  Since  ~2Jf£ 

the  argument  of  the  last  for. h in  (4.27)  lies  in  the  4th  quadrant 
with  phase  between  ~ and  O . 'fhe  ionh  of  this  argument  then  has 
the  phase  limits  to  .028  radians  (1,6°).  This  can  be  shown, 

for  example,  by  the  general  relation 


sZ±J*2±iZ_  fg 

Ctift  2 t j.  coj  2 y 


(4  Z 3) 


Multiplying  by  Zo  with  phase  i77/)) , then  gived  Z the  phase  O -to  .6/3 

radians  ( 46.6°)  (The  upper  value  of  this  phase  comes  when 

('in  2 hi/imr.  2 b t ) which  gives  the  phase  of  tori/>(-j9tJl) 

has  its  minimum  value  at  2 A Z ® 3.  S3  Just  less  than  •£pr ). 

For  the  case  )io)  s’  / 


fonh  x ' j 0^x2^ 

(the  sign  of  the  imaginary  term  depends  on  the  breach,  but  does  not 
matter  hereL  Then  (4.27)  becomes 


Z = 2*  [x  + ki  -*(&£+  @)J  - 2a 


S/S)h(2*f2k£)+j  s)n  2t?Z 
ccshUxtZ&J)  - cos  Zkl 


(4Z9) 


When  4'*^,  corresponding  to  £'/^0  infinite,  and  /t£  is  small, 

(429)  has  the  form  2?  * Z»  (f+h/fti,  while 
for  forge  2)  , Z *Z0  Thus  Z ranges  >n 
phase  from  f71/^)  to  .758  radians  (4-3.4'). 


Hence  the  transformation  of  any  Z ' with  phase  (”/*)  , down 
any  length  of  line,  results  always  in  an  impedance  i-em&ining  in  the 
first  quadrant  and,  in  fact,  rapidly  approaching  the  value  2>  of  the 
transforming  line.  This  is  shown  explicitly  by  the  diagram,  where 
impedance  spirals  starting  from  various  values  in  the  first  quadrant 
are  shown  rapidly  converging  on  Z , and,  of  course,  always  remaining 
in  the  first  quadrant. 

This  result  suggests  the  general  theorem  that  all  thermal  impedances 
seen  by  a thermal  wave  propagating  in  a rod  of  arbitrary  cross  section 
and  material  with  an  arbitrary  termination,  lie  in  the  first  quadrant, 
hence  are  resistive  and  capacitative . low  any  impedance  in  the  first 
quadrant  may  be  realized  by  proper  choice  of  Ho  'y  H©  and  -& £.  , as 
shown  by  the  discussion  above  of  {4-  c.T  ) and,  particularly,  by  the 
diagrams  of  the  impedance  spirals  produced  by  transformation  down  a line. 
Hence,  with  this  theorem,  we  have  a complete  existence  theory  for  thermal 
impedances , 

Such  a result  is  nr.de  plausible  by  the  argument  that  any  termination 
may  be  considered  approximately  as  a series  of  thin  slices  of  uniform 
cross  section,  hence  as  a series  of  sections  of  thermal  transmission  lines. 
Evidently,  if  we  start  from  some  boundary  condition  such  as  7 ~ ~ Oj 
or  ^T/$r)  -O  , corresponding  to  impedance  & ar 
the  result  of  transforming  through  tlie  series  of  transmission  line  sections 
remains  in  the  first  quadrant  by  the  above  discussion.  Even  a radiation 
boundary  condition  (into  a vacuum)  in  which  ~K3T/sn)  - /V7” 
where  H is  a positive  constant  . II.  larsiv..  and  J.  Jaeger  - Conduction 
of  Heat  in  Solids,  Oxford  1947  p«13)»  ui ay  be  included,  since  this  corresponds 
to  an  impedance  of  (1/H)  (per  unit  area)  which  also  starts  in  the  first 
quadrant . 


This  argument  is  difficult  to  aake  rigorous  by  passing  to  the 
limit  of  zero  thickness  slices  because  it  requires  consideration  of 
the  discontinuity  impedances  introduced  by  the  successive  changes  of 
cross  section  going  freo  slice  to  3 lice  . These  impedances,  which 
are,  in  fact,  just  lumped  series  impedances,  as  shewn  in  section  5* 
do  go  to  zero,  as  the  slices  get  thinner.  To  neglect  them,  however, 
they  really  must  vanish  to  higher  order  than  the  thickness  of  the 
slices.  This  does  appear  reasonable  from  the  equations  of  section  5 
for  the  change  of  circular  cross-section,  where (5. 33)  for  Z.  and 
(6.6)  for  CcJ  show  that  as  f , the  ratio  of  radii,  approaches  / , 

the  factors  ) appearing  in  each  term  of  iv.  approach  zero 

as  (T-/)1  j since  U,  (j8n)-o  . However  is  only  the  first  term 

of  t he  variational  expression  for  H , and  consideration  of  higher 
terms  would  be  more  difficult.  In  addition,  even  if  uhis  can  be 
carried  through,  actual  abrupt  changes  of  cross  section  could  still 
be  present  in  the  termination  and  the  theorem  would  have  to  be  proved 
separately  for  them.  Accordingly  we  look  for  another  type  of  proof. 

(4g)  General  Pr  of  that  all  Thermal  Impedances  lie  in  the  first  Quadrant  „ 

First,  we  note  that  the  restriction  of  the  impedance  implies  a 
restriction  on  the  reflection  coefficient  & of  a wave  incident  on 
that  impedance.  If  Sr  is  the  terminating  impedance  and  So  the 
characteristic  impedance  of  the  line.  (4.19)  gives 

/ + & £ _ zrg  - So  ( *rv  3 OJ 

/—  f2  ' + So 


The  transformation  (4*30)  takes  the  complex  impedance  plane  conformal' 
into  the  complex^  plane,  with  the  first  St  quadrant  going  into  a 
portion  of  the  unit  circle.  Thus  putting  7*0  - /■+*.  , which  loses 

no  generality  but  just  fixes  the  scale  of  impedances . 


(4®30)  becomes 

2 - Ll  ~ (e-rtf 

(>-ieis‘e)'~  1+  K!1  -£-«• 


(4. 30 


where  the  conjugate  oi‘>_  in  denoted  by  <2*  ~ j£je  * , In  (4=31)  the 
denominator 

ft! Oil1-  ft-  ft'  * //  //?/"“  1 1 Cos dP  - (i~  /til Cos &)Z*  {£?S;nl6>  *0  (*■  32) 


.««(>■ -T^V 
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',ence  is  always  positive,  Then  if  Z.  lies  -211  the  positive  real 
axis,  we  must  have 

/ - /^y *- 4 (fi-  ft *)  - /-  /*?/*  *■  z /n * s°  (+■  ~3<*J 

j 

j 

/-  /£/<  +j(#  -&*)  - / - /#/Z  " -?  JtfjS/nQ^o  (4:33b) 


Putting  (*  1J<7  } in  {4  33b  ) gives  - 4/#/sss>0^o , hence  j \ 

an<4  solving  (4.32*)  for  /<£/,  t 


/&  / z V f'n  i & + S'n& 


-77  < & <-  O 


(4.3*) 


(4«34)  gives  the  locus  oft-  when  Z±  is  real  and  positive,  and  evidently 
lies  on  an  arc  dthin  the  unit  circle  stretching  from  4/  fo  “/ 
below  the  real  axis 

Similarly,  when  Zb  is  positive  imaginary,  the  conditions  (4*  33cj) 
o/)j£>BX9  interchanged,  and 

/<??/=  ffrs'in*  O > Oi&±rr  (4-5S) 

which  is  an  arc  obtained  fcy  reflecting  ( 4 »34)  in  the  real  axiso 

Evidently  the  condition  that  «?*  lies  in  the  first  quadrant,  requires 
inequalities  fear  both  ( 4.33  ) a and  b,  and  leads  to 


f c )/  y-  s// 7 1 (9  ~~  & <0  £ & 4 7t~ 

I if/r  sin  ‘ aT  t*-  sins  -rrceiO 


(4.36) 


so  that  n.  lies  in  the  elliptical  region  in  the  unit  circle  between  the 
arcs 

(4.35)  and  (4.36)  . The  transformation  is  illustrated  in  the  diagram 

, . L Cfoti) 
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6(0,  fZ  i) 
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4 
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£(<+  fZ  ~U) 
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Thus  the  theorem  will  follow  if  we  can  show  generally  that  the 
temperature  fielu  in  the  rod  bias  a principal  mode  of  the  form 
r ' +-  ? s-~l *■'*' , (talcing  the  termination 

at  b-o)  where  is  restricted  in  the  above  fashion. 


In  the  case  of  electromagnetic  wavesin  waveguide  or  acc-^tical* 
waves  in  pipes  a restrict?. on  on  \ is  obtained  immediately  fr cm  con- 
servation of  energy.  Thus  since  the  reflected  wave  cannot  carry  away 
more  energy  than  the  incident  wave  carries  into  the  termination,  and 
since  the  energy  flow  is  proportional  tc  the  absolute  square  of  the 
amplitude  of  a prop-gating  wavs,  this  gives  immediately  /..€/ - £ /, 
and  R lies  within  the  unit  circle,  or,  correspondingly,  Z,  lies  in 
the  right  half  piano  ’Then  is  real.  For  thermal  waves,  this 

argument  fails,  since  thermal  waves  do  not  cany  any  energy  ( the  time 
average  heat,  flux  is  zero)  and  would  in  any  case  be  inadequate  since 
« is  restricted  more  than  by  /<*r/  ^ / , All  that  one  can  say  by 

direct  physical  argument  is  that  the  thermal  wave  represents  propagation 
of  a disturbance,  and  we  expect  that  a disturbance,  which  here  propv^ates 
by  a diffusion  process,  must  keep  diminishing  as  it  spreads  and  will  be 
reflected  in  general,  with  diminished  amplitude. 

However,  the  formal  proof  of  the  onergy  theorem  from  the  wave 
equation  and  boundary  conditions  is  obtained  by  applying  Greer's 
theorem  (4.6)  to  the  tangential  electric  field  and  its  conjugate.  We 
try  a similar  procedure  here  applied  to  T and  T \ Noting 
that  T - -rt  lT,  C (4.6)  gives 


(jtL-  *')( IT/1 4 1/  ' -/rr*F  - T’  ) <*# 

'JJ  A-  J > ds- 


(4  37) 


where  the  volume  V includes  a section  of  a rod  leading  into  an 
arbitrary  termination o Cn  all  boundaries,  either  T~o  »/•  tJTln) 
except  on  the  cross  section  .*  of  the  feeding  rod,  hence  only  that 
surface  integral  remains  in  (4.37) > and  there  %r>  - - . 


SMS  UUW  pUO 
. if  3.  ^ 

Jz.fr  * r Si  & 


" h Z e Jo  ( 


(4.38) 


as  a general,  exact  expression  for  / within  the  rod  in  which  only 
a principalmode  is  incident  giving  rise  to  reflected  principal  mode  and 
higher  mode  waves.  The  reflection  coefficient  is  r at  cros3  section 
o z o > which  could  be  at  tiny  position  along  the  rod.  Then  the 


termination  is  defined  ao  the  materlalfor  £ Po  , and  ■)'  is;  the  reflection 
due  to  that  termination,.  Although  for  convenience  we  have  used  circular 
cross  section  higher-  modes  (with  cylindrical  syirmo try)  in  ( 4,36  ) this 
is  not  an  essential  restriction  since  the  transverse  parts  of  the  solution 
of  the  wave  equation  in  any  cress  section  fora  a complete  orthogonal  set  of 
functions  in  that  cross  section,  with  discrete  real  characteristic  v?  me?  <C? 
which  appear  in  the  ->?/»*/» </<>,-7<r  through  Pn 1 - J**-  *7-> \ 

Returning  to  (4.37),  we  put  X - ( f-t  / ) /■  and  breakup  the  integral, 

over  IX  into  an  integral  over  the  termination  and  one  between  3-  crs>t/&, 

to  obtain 

4Ax[flr/lxix 

7e-r/s<oioti-jii 


+ faf  XA/T/f  : -JX4- 

3 > J * 


£4 


:riL 


7~- 


(4.3?) 


New  from  (4.38) 

pA/T/l  = 1 V at e'****  #*ezi**+A<Lt /Anfe 


1 ~<( Xn~£n)j-  . 

* T*~  . 

(4  4o) 


where  use  has  been  made  of  the  orthogonality  properties  (5, >7)  of  the  mode 
functions  in  the  cross  section,  and  of  ^ 

Al3c  from  (4.38) 


— , , -J  (#„-*/)  3.  - 

y sJ^.  ■ nt  1 

A 

(44/) 


Hence  (4»39)  becomes,  after  integrating  (4.40  ) frem  j?- 
noting  that  the  principal  mods  terras  cancel  out, 

/-/4i/ 1 -M#  -#  *)  - 3±  f/T/~'ji/  / y y 3.. 

4 / -a(j<a  -43) 

/h-.-.-n  matron 


to  O’ 


and 


(4  42) 


( 4 c. 42)  is  independent  of  ^ and  also  follows  by  putting  in  (4,39) , ( fcoAl)  , 

r-r  xn  - "On"] ^ where  the  root  wit  h positive 


is  chosen  ( so  that  e"  ~n*  ~AU> 
+ ).  hence  X0 

small  for  large  n )„ 


imaginary  part 

! t-tcnr.  -?ci  in  the  direction  of  propagation 
must  also  have  a positive  real  part  ( which  is 


Then 


(j£n  + X'4)P’0  ; (<£"  ~ Xnr  ) XD 


(4.43) 


and  (4.42)  leads  to  the  conclusion 


/_  /£/  v (£  - & *)  s / - /&/*  * Z ^ d? 


(?A+J 


Comparing  ('vo44  ) with  (4»31  Jshcws  that  we  have  proved  that 

always  has  positive  imaginary  part,  or  lies  in  the  right  half  plane. 

For  -rr&  s .4  O , this  also  gives  At  Zt  since  /te  £*  ■ /-/*/ 

~4/vi/ Z/P&  on  using  {4- 49) 

However  for  o<etrrt  /?«■  Jt.f  could  bo  negative;  for  example  if 

8-n/LJ/fi}~fZ+l,  then  (4.  ft)  /J  so+Ji-fivd , and  /??  £?  - — df  / r/T) . 

Thus  only  part  of  the  theorem  restricting  Z?  to  tho  first  quadrant 
has  been  proved. 

The  other  condition  or,  F*  is  obtained  similarly  by  applying 
Green's  theoremto  (°%>})  ond  Ik  place  of  (4.3?)  we  obtain 

(* •••  yw * *m&W)  - W-  a $0 

(4  4S-J 

whore  the  reduction  of  S to  just  the  cross-section  V depends  on  the 
vanishing  of  H>n  ^ ) when  either  or  {JT/m)  -o . 

Thus  vdien  7=0;  (derivatives  in  all  tangential 

directions ) , hence {s"r4t^) * - se * 7“  - £*•(* -)  - a. 

Therefore  cos(f>.^)('J'r/»rd)  = ot 

where  <e;(o.7)  is  the  cosine  cf  the  angle  between  the  > axis  and  rr-.<- 
normal.  In  the  other  ca3a  (sVdn)  -o  , hence  (? « Then 


) _ ZL(-ZZ-~  r o ' ( tan,  2 j ) - 7Z.Z1 — cesff /’/>  //,'/ 
dnidtern  ' ~ dnJfvn 


In  : lace  of  (4.39)  we  now  hav© 


Sitf/gftr  +&>£*/#/*}* 


£T  sdlZ'*) 
> / 


and  putting  (4.33)  for  7”  in  (4.46)  gives 

/-  /e/'+.i(e-K'j  = 4 *i  tk miZMB) 

A 7.  T°rn-..r«,n°s>  * /?;/ 


^*4 
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Then  (4o43)  and  (4.47)  lead  to  \ 

/- /ft)'4-*) (£- £ ~)  s /-/^/c  ~ 2/#/s;, 7 


(+  4&) 


which  with  (4*31)  and  (4.44)  completes  the  proof  that 
lies  in  the  first  quadrant. 


always 


The  significance  of  the  conditions  (4.4S)  and  (4.44)  is  further 
brought  out  by  the  behavior  of  the  magnitudes  of  7T>  and  )/>  as 

functions  of  j . 77  and(*%j  )p  are  the  lowest  mode  parts  of  77 

and  are  also  the  mean  values  of  7~  and  over  the  cross  section 

in  view  of  the  orthogonality  relations  of  the  higher  mode  functions. 
(!>7%j.)tc  is  thus  proportional  to  the  net  heat  flow  through  a cross  section, 
as  shown  in  (5.14).  <x, 

It  follows  easily  from  -a  <6  *•  that 


•1ST ^ -zk/**m*-f*r*z**+i  to* 

/> 


'*‘*7  _ 


(4.  19) 


£ y^/%  (ft  V“*)J  (4.^0) 


Thus  at  £ 'o J (4.4#)  and  (4.44)  reduce  to  the  conditions  that 


{*.  *y? 


Since  ) was  an  arbitrary  position  in  the  rod,  this  means  that  the 
mean  values  of  77  and  ( J )*>  continually  decrease  in  magnitude 
as  wo  move  away  from  the  source,  and  nothing  like  a node  or  minimum 
in  the  magnitude  of  the  field  can  occui  - Thus  a resonant  condition 
and  true  standing  waves  can  never  occur  for  temperature  waves,  as  we 
might  also  expect  from  the  fact  that  only  ^apacitative,  but  no  induc- 
tive reactances  are  available  for  these  waves.  Hence  no  resonant 
combination  with  aero  reactance  can  be  put  together. 

A direct  proof  that  (4.51)  holds  may  also  be  developed  by 
relating  the  values  of  77  and  ( 3Tj0/jj-)  in  the  rod  to  those 
at  the  boundary  of  tne  tcnuizuition.  One  considers  the  mean  values  of 
T and  (*7/^)  over  croos  ?=;octions  of  an  extension  of  the  rod 
into  t^e  material  of  the  terKiiisvcion.  A theorem  is  easily  proved  that 
shows  -nese  man  value 3 can  never  have  a maximum  in  the  conducting 
material.  Hence  if  the  viluos  increase  going  into  the  material  at  the 
bc;uJ,U'>,  'Ir  ".-.l  inoresse  ail  the  way  to  the  source,  which  gives 

b»aGdaiy  if  start,  for  e-xample,  with  th©  condition 
= J th&a  /T/  r.-.r.t  increase  on  going  into  the  material. 


i 
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vhe  theorem  about  the  maximum  follows  simply  from  the  identity 


' 0>  ‘cTTlTJu 


(4-T2) 


where  7~f})  is  the  mean  value  of  T(j-)  over  the  cress  section  of  the 
rod  or-  the  extension  into  the  termination,  (4.52)  comes  from  Green's 
theorem  applied  to  77; ) and  1 ( the  cross  sections  must  lie 

completely  in  the  material.) 

Various  complications  derelop  when  the  cross  section  intersects 
the  boundary,  and  in  the  consideration  of  ) with 

boundary  condition  7*  ° and  vice  versa,  that  make  this  argument  difficult 
to  carry  through,  and  the  method  given  above  seems  preferable. 


5.  THE  CHANG-:-:  CF  CRQ.S  SECTION  FOR  CIJiCULAP.  ■ CYLINDE.TS- 
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General  Outline  and  Background 

W«*  now  solve  explicitly  the  rrobleia  of  the  temperature  distribution  v round 
a June Lion  between  circular  and  cylindrical  rods  of  different  diameters  and  materials. 

As  indicated  previously  this  will  separe  a into  a static  and  a dynamic  problem®  We 
are  primarily  interested  in  the  dynamic  problem,  namely  reflection  and  transmission  of 
a principal  mode  temperature  wave  at  the  junction,  since  it,  is  this  wave  on  which 
measurements  are  made  in  the  experimental  technique  for  determining  diffusivities. 

In  the  following  we  separate  the  static  and  dynamic  problems.  The  latter  is 
solved  by  first  expressing  the  temperac-ure  fields  in  characteristic  mode  expansions, 
then  matching  the  expansions  at  the  junction  plane  z»0  and  obtaining  an  integral  equation 
for  iii  the  cross  section  of  that  plane®  Voltages,  currents  ana  characteristic 

impedances  for  the  principal  mode  behaviors  in  the  two  rods  are  defined  in  accord  with 
tire  transmission  line  analogy  of  the  last  section.  The  integral  equation  leads  directly 
to  a variational  expression  for  a series  impedance,  j?  , representing  the  effect  of  the 
junction  plane  on  the  principal  mode  amplitudes,  i.e.  voltages  and  currents,  in  the  rods. 
Expanding  the  function  ( ) in  modes  of  the  first  rod,  substituting  the  expansion 
into  the  expression  for  <?  and  usung  its  variational  character,  gives  an  infinite 
set  of  linear  equations  for  the  mode  amplitudes.  In  turn  2 is  expressed  in  terms  of 
these  mode  amplitudes.  The  numerical  evaluation  of  2 , for  appropriate  constants, 
is  then  carried  out  in  the  next  section. 

Analogous  electromagnetic  problems  are  concerned  with  propagation  of  radiation 
in  a wave  guide  and  through  a wave-guide  junction.  The  formulation  of  these  problems, 
which  introduces  a transmission  line  analogy  and  obtains  variational  expressions  for 
elements  of  the  equivalent  circuit  representing  the  junctions,  is  due  to  J . Schwinger 
Later  work  along  these  lines  is  also  relevant. 


(2)  J.  Schwinger  - Unpublished  lectures  at  the  K.  I.  T.  Radiation  Laboratory  , 1943°*46 „ 

(3)  D.  Saxon  - "Notes  on  Lectures  by  J.  Schwinger,  Discontinuities  in  Waveguides", 

Feb.  1945  (develops  transmission  line  and  equivalent  circuit  analogy,  solves  a 
number  of  simple  discontinuity  and  junction  problems  by  variational  formulation 
and  the  static  method,  including  capacitive  diaphragms  and  the  change  in  height 
of  a rectangular  guide.) 

(4)  "the  Waveguide  Handbook,’  edited  by  N..  r".rcuvit.;?,  McGraw-Hill  Book  Co.  (1949) 

(tabulates  equivalent  circuit  elements  for  a great  variety  of  waveguide  discon- 
tinuities and  .junctions,  discusses  the  transmission  line  analogy,  but  does  not 
give  any  mathematical  methods.) 


tV 

■j  (5)  L®  Lewin  - Advanced  Theory  of  Waveguides,  Lliffe  and  Sons,  Ltd.  London  (1951) 

(Presents  some  of  the  mathematical  methods  of  solving  some  of  the  discontinuity 
problems,  including  E plane  and  H plane  steps;  gives  expensive  bibliography) 

4 

(6)  J.  Miles  - J.  Accost.  See.  .bn,  16,  1*,  (1944);  17.  279  (1946)  (sets  up  the  problem 
of  reflect iou  of  principal  mode  sound  waves  in  a circular  tube  at  a change  of 
cross  section,  giving  variational  formula  for  equivalent  circuit  element;  no 
numerical  results  or  calculation  nethods  given) 
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The  present  problem  di  fere  from  the  electromagnetic  waveguide  problem 
( and  the  acoustic®!  problem  of  sound  propagation  in  pipes)  in  the  following  ways. 
First,  the  quantity  sought  is  a scalar  field  ( as  in  the  acoustical  case)  rather 
than  a vector  field  ( in  fact  the  electromagnetic  problem  seeks  two  vector  fields 
E and  H)  so  that  the  differencial  equations*  characteristic  mod os  and  boundary  con- 
ditions are  *11  much  simpler.  (Ultimately,  however,  only  problems  reducing  to  one 
dimensional  problems  in  a junction  plane  or  on  a surface  of  discontinuity  cen  be 
solved,  in  the  electromagnetic  case  as  well.)  Second,  the  principal  mode  field  is 
constant  over  the  cross  section  and  has  no  cutoff  (l.e.  minimum)  frequency,  a situ- 
ation analogous  to  the  acoustical  problem,  but  unlike  the  electromagnetic  case. 

Third,  all  modes,  including  the  principal  mode,  have  complex  propagation  constants, 
so  that  even  the  principal  mode  attenuates  and  the  corresponding  transmission  line 
has  a complex  characteristic  impedance.  This  situation  is  unlike  both  the  electro- 
magnetic and  acoustical  problems.  Fourth,  the  actual  physical  problem  here  separates 
into  the  superposition  of  a static  and  s.  dynamic  problem,  w tn  solutions  of  quite 
different  character,  of  which  only  the  la  \,er  is  represented  by  the  transmission 
line  analogy.  Finally,  we  note  that  the  problem  of  the  change  of  cross  section,  as 
solved  here  by  th6  use  i f a variational  expiree sion  for  the  equivalent  series  impedance, 
includes  the  possibility  of  on  abrupt  change  of  material  at  the  junction  plane,  a 
generalization  of  more  importance  for  this  problem  than  for  the  electromagnetic  or 
acoustical  problem. 

(5b)  Description  of  Physical  Situation  and  Form  of  Solution 


As  in  the  diagram,  rod  I with  radius  xl  and  constants  it  and  D is  at 
the  left;  rod  II  is  at  the  right  with  the  corresponding  quantities  R ' f it’  ,Q' 
designated  by  primes  as  in  Section  3?  rod  II  will  by  convention  always  have  the 
larger  radius.  The  junction  pLane  is  at  z , where  the  changes  of  cross- 

section  and  material  take  place. 

To  illustrate  a definite  physical  situation  we  assume  a time  varying  source 
of  heat  at  a free  snd  at  - L,  in  rod  I,  sufficiently  far  from  the  junction  so 
that  only  the  principal  mode  cornniunica  os  between  source  and  junction,  and  let  rod  II 
extend  to  infinity.  Then,  as  shown  in  Sect  on  2,  the  solution  to  this  problem  is  the 
superposition  of  two  problems.  One  is  a static  problem  (actually  quasi-static)  with 
a steady  heat  source  <?.  ot £ ~ and  a steady  heat  current  flowing  into  and  past  the 
junction,  while  T satisfies  Laplace's  equation.  The  other  is  a dynamic  problem  with 
source  at  / - - 4,  , and  T satisfies  the  wave  equation  (V*  *-rtl ) 7~ * O, 

Restrict  Inge  a*  attention  for  now  only  to  the  dynamic  problem,we  expect  T 
in  rod  I to  consist  of  travelling  principal  mode  waves  in  both  directions,  and  in 
rod  II  of  outgoing  principal  mode  t/aves  only.  The  situation  in  rod  I is  described 
completely  by  the  reflection  coefficient  & seen  by  a wave  approaching  the 
junction  (when  rod  II  has  a semi-infinite  or  matched  termination).  /?.  can  be 
calculated  from  a knowledge  of  the  equivalent  circuit  of  the  junction,  as  shewn  in 
Section  4°  The  f olds  in  rod  I then  take  the  Tones 
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which  agree  with  (4.4),  have  voltage  reflection  coefficient  /? 
and  satisfy 

r^tfiZ)  - -KAU'Q'iXZ*  Q< 


at  ^ -< 


rx-2; 


fc2i*LZ<</,  fu  f)  takes  The  form 


ixc, 


Vis)  * 
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XAX2, 


fry 
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(5«3)  is  i(.;; Mediate  if  the  assumption  is  made  that  the  reflection  at  the  source  end 
of  the  reflected  wave  can  be  neglected „ 

We  proceed  now  to  the  rigorous  formulation  and  solution  of  the  problem, 
which  will  lead  to  a circuit  representation,  hence  to  the  reflection  coefficient  of 
the  Junction  0 

The  Mode  Expansions 

Expand  (^TOj.)  in  rods  I and  II  in  the  complete  3et  of  modes  with  circular 
synmetry,  as  found  in  (2.12),  but  taking  only  principal  modes  incident  from  left  and 
right . Then 


(jff  - *S>n)tS  ' 1 A',  J,  frfr)?  iXn> 


fry 


n *o 
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where  we  have  introduced  only  higher  modes  which  go  away  from  the  junction,  and 
jCj  ff, , are  defined  in  (2.9).  (Z  tt).  ,e 

note  that  {dT/^)u  on  the  annular  rnetal  surface  between  R and  R*  in  the 

plane  z«0»  Hence  in  that  plane,  values  of  exist  only  in  the  circular 

region  (radiu  R)common  to  the  two  rods.  This  is  why  the  basic  expansion  (5.4)  is 
given  for  (*T/>  s)  rather  than  for  T.  New  the  values  of  at  the  junction 

satisfy 


(s.s) 


since  heat  flow  across  the  junction  is  conserved,  Cail  the  function  in  (5.5) 
E(r) . Then 


£(r)  ~ A' 2.  Ai  J;  V 

/y  -o 


/C2  Jo(C,’rJ 

r>:o 


(S'.  6) 


Now  making  use  of  the  orthogonality  properties  of  the  (V*  fj  f namrly 

J Zi (£> r)Ufa r)r</r  r f ?±r U frT.  jL U?  rJ ~ -- j*o  for  rr,*n 

* /rn  " ' J"r‘~ 
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we  can  evaluate  /tp  onct  by  multiplying  (5.6)  by  r-A(fn  r)  or  roc-  (fr  r) 

and  integrating  from  0 to  R or  0 to  R1  respectively,  giving 


A >wC^V^'/T  M WU I dr  J A"  = /r (£/?)/ ^fri  J°  (r» >>ra'r  (S,&) 


Integrating  (5.4) 
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where  the  constant  of  integration,  independent  of  is  discarded  because  it 
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does  not  satisfy  the  differential  equation  ( ( 7 - <7>J, 

Substituting  the  lowest  mode  terms  in  (5.9 ) by  VO>/c  , ^ fj-Vc  <}'*■*& 

^ r-  J42/  VI*.).  - —Cos  Xs..  ~ s4,€>~A*J 

C ' c.  • * ^ TJi 
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and  similarly  i**  (5°4)  introducing  / V j)  t / 'fj ) g/res- 
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(5d)  Fixing  the  Characteric  Impedance 

The  constants  Cj  Z„ , £T  ^ ZJ  are  thus  far  arbitrary  in  the  formulation 
of  the  transmission  .line  analc  £j.gc  for  X*  ods  1 ana  II.  Simple  additional  requirements 
will  fix  their  values,  however.  First,  vs  wish  to  represent  the  effect  of  the  dis- 
continuity by  a lumped  circuit  inserted  between  the  two  transmission  lines.  However, 
a lumped  circuit  necessarily  implies  that  the  linear  relations  between  V,  VJ,  /,.  TJ 
satisfy  reciprocity,  hence,  as  shown  in  fee Lion  (4b),  we  must  have  (4.13), 

K/l  it  Za  _ SfA'  Zf 
= ' (C  )L 


(s:  /2) 


Second,  we  wish  to  simplify  the  equivalent  circuit  ar  much  as  possible. 
Now  we  see  from  (5.11)  that  by  choosing  the  natural  reference  plane  J ~ 0 at 
which  to  evaluate  currents  and  voltage.  \>v,  can  make  J0  * J,  ‘ for  all  conditions 
of  excitation  (any  combination  of  principal  mode  waves  incident  on  the  junction), 
since  at  ^ an 7 ■*c  drop  out  of  the  equations  for  I(j)  and  I*  (a;.  Then, 

taking  /?.  and  A*  from  (5.S)  and  noting  that  the  common  factor,  jo  ’'<Efnjrc/r 
cancels  out,  we  obtain  from  the  requix-emant  Z,  - 7. * 


/r/?jg  2> 

~~c~ 


/<  >?'  X4C'&o 
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(r-  n) 


Now  A --77r<ly  A ' - 77 £? 


. hence  dividing  (5.13)  by  (5.12)  gives  C-  - C ' y 


a 


/ 

.id  thus  J 
particular  rod. 


a universal  constant,  independent  of  ulte  a -sciai  properties  of  a 


(5  .12)  and  (5<>13)  are  then  satisfied  by  choosing  St  proportional  to 

( )■  The  proportionality  constant  may  now  be  fixed  by 

the  requirement  that  I(  ) be  simply  .hr  total  hea  . flux,  <2  , passing  through 

any  0/.-OS3  section  of  the  rod*,  We  liave 


- /ZA-2  2?  T 

'"lz 
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( r.  /*) 


on  using  (5.H)  and  noting  that  by  (5.7 ' / ^ so  that  only 

the  principal  mode  contributes  to  the  ho/i^  flux  through  a cross  section0  Hence, 
if  we  choose 


2 

STAk 


c --/ 


we  have  the  simple  relations  from  (4*3)  and  (3*14) 

I //>;*  T,(j) 

id)  -■  -K'+af),  = 


(s:  /<s) 


where  7#  ((3)  and  ( ^j)/?  refer  to  t!ie  principal  mod?,  parts  of  7f£)  bxa(^ 

Since  Q.  must  be  continuous  across  each  cross  section,  including  the 
one  at£a^  s (5»14)  and  the  requirement  that  I Co)  be  continuous  lead  immediately  to 
(5*13)  and,  in  fact,  provide  an  alternative  derivation  of  (5  A3)® 

The  continuity  of  2 /'o)  means  that  the  equivalent  lumped  circuit  is  simply  a series 
impedance,  with  no  elements  in  shunt  across  the  line,  and  explicit  expressions  for  the 
series  element  will  be  derived  below  B Note  that  it  is  not  possible  to  make  I /foj 
continuous,  since  by  (5  AO)  the  coefficients  oc  and  =c  ' do  not  drop  out  at  £ = 
and  these  vary  with  the  excitation  conditions,*  Honour  7~  C °)  must  be  continuous  over 
the  region  o C r < /?  y sc  that  we  must  have  , by  (5  AO) 


~rJ  (o)  = 7“  71  Co) 
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(5e)  The  Equivalent  Series  Clement  and  the  Variational  Expression 


Now  the  discontinuity  in  /-«*  - Yo‘  , must  bo  proportional 
to  the  single  value  of  Mo)  which  measures  the  lowest  mode  amplitude  of  E(r),  since 
all  equations  of  the  problem  are  linear.  Increasing  E(r)  by  a factor  will  increase 
I/O) and  Vfo)  -\/‘fc)  by  the  same  factor.  Hence  we  may  set 


1 yfoj  - \Z‘(oJ  « 2 Ifo)  --  - 2 rrZj^rS  fr)  or 


(S-./3) 


thus  defining  an  impedance  2 which  is  independent  of  the  excitation  conditions  and 
determined  the  overall  effect  of  the  obstacle  on  the  lowest  mode.  Then  (5*17)  becomes 


/ f 2 c~  ^ 

* Zrrf£(r)r  Or  / /T/?1 

ypz  ^ ^ -6 ; rv* ' /?  :/ 


4 (fr>r)rOr 


(5.19)  is  a homogeneous  integral  equation  for  E(r),  for  which  is  the  eigenvalue. 

The  relation  (5.13)  clearly  has  the  following  circuit  representation 


A 7.' 


This  has  a single  series  impedance  of 
magnitude  2 in  the  line,  as  ex- 
pectec  from  the  continuity  of  13. 


From  (5.19)  we  obtain  a variational  expression  (i.e.  stationary  for  snail 
changes  of  E(r)  ) for  2 in  terms  of  E(r),  by  multiplying  by  & tr)rcfr  and 

integrating  over  o to  /?y 

Z r 7 -/ ~ , l—  2.  ‘ ■ (fsfoJo  (fr.  rjrc/r ) 2 Zj/T/ ; (27#  J (fsf'M (C*  r)retr)  j 

c n/JtrtrJr] ‘■[K /?•$■.,  / /•/?'*,,.  .r.u  /J 

(S'.  20) 


The  statior.ary  character  of  2 is  easily  proved  by  "differentiating"  (5.20) 
with  respect  to  £(*")  (i.^.taking  the  functional  variation  of  £(r)  )c 

Thus  we  write  (5*19)  in  the  schematized  form,  (the  stand  for  the  constant 
coefficients  in  (5.19)) 


'* 


and  (5.20)  in  the  form 


(f*£(r)rctr)1  , . 
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Taking  the  variation  *E 


of 


of  P 
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(5.22)  becomes 
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£:  and  <f£ 

(f  Z £ (r)  rc/f)  “ + 2 2 ( fz  rats t )[  f-S£  rcfr)  *2  <i  /%  ( jtU*  rdr )(fr£<Je  r Jr) 

which  gives 

J2  ( fZ rdr)1  - 2 fj£rc/r £ j?  f£  rdr  y-  £ Jc  f ^Uo  rdr]  -o  (S'.  2 3j 


where  the  right  hand  side  of  (5*23)  vanishes  by  the  original  integral  equation  (5.21) 


We  now  rewrite  (5.20)  in  terms  of  the  coefficients  of  the  higher  modes 
using  (5.8)  and  (5.6).  We  have 
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hone©,  using  (5.7) 
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and  putting  - f 

fixes  the  magnitude  of  £(r)  ) 


(which  implies  no  loss  of  generality,  but  simply 
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Putting  (5*30)  back  into  (5.29)gives 
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• ( 5f ) Linear  Equations  for  Mode  Amplitudes  and  Jimplifxed  ;-  orm  for  Series  ^leggnt 

Differentiate  (5.29)  with  respect  to  the  sy>*t2,~  '■*>  to  give  a set 

of  linear  equations  for  the  Am 
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Trie  equations (5. 30 )are  an  infinite  set  of  linear  equations  for  the  coefficients 
of  all  the  higher  modes  of  E(r ),  which  when  solved  and  substituted  in  the  variational 
fora  for  5?  , give  by  (5»32)  a series  for  The  first  approximation 

for  J?  is  simply j?»  Coo/krr/fi  and  the  second  is  obtained  by  solving  for  >?, 
in  (5*30)  with  <4,  />;T»o  , giving 


/C 77 


fc. 
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mm<L 


(5-  33) 


This  ccoipletes  the  formal  solution  of  the  problem  and  we  are  left  with 
Just  the  problem  of  numerical  evaluation  of  the  Cm/t%  which  are  then  substituted 
in  (5.39)  and  (5.32)  or  (5.33).  In  Section  6 the  expressions  are  repeated  in 
appropris.ee  form  and  numerical  evaluation  made  » 


We  note  that  after  2F  is  found,  the  reflection  coefficient  is  ob- 

tained from  transmission  line  theory  as  expire  sed  in  (4.20)  - (4.23),  but  now  we 
deal  with  the  simple  case  - &> , 2,,  — 2 2 j?/4.  * 2 
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where  i?  4 is  the  termination  of  line  II,  as  viewed  from  the  Junction . 


Mode  Amplitude  Equations  by  Direct  Matching 

It  Is  worth  noting  that  the  equations  (5.30)  for  the  /h*  obtained  from  the 
condition  that  B be  stationary,  also  follow  directly  from  the  temperature  matching 
equation 


To  eee  this  we  expand  the  integrals  of  E\'rJ  us/n§  CST6J 
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Hence  (5.19)  becomes 
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Multiply  (5-35)  by  J. Y J t and  integrate  from  o tv  ft  to  give 
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which  agrees  with  (5»30)o 

Thus  w* > check  the  algebra  leading  to  (5^30)  and  prove  again  that  /L  is  stationary. 


6,  v mJKF.RICAL  EVALUATION  OF  THE  SERIES  ELEMENT. 


AND  REFLECTION  COEFFICIENT  &. 


(6a)  Summary  of  Procedures 

The  computation  of  the  equivalent  series  element  2 for  the  Junction  of  the  t.vjn 
rods  is  still  a complicated  matter,  since  the  Bessel  function  summations,  which  give 
the  Cgm  must  first  be  evaluated*  These  sums  all  converge  for  large  n like  ((4>)*  i 
but  may  be  slower  for  a practicable  number  of  terns*  Accordingly  the  convergence  is 
assisted  by  summation  formulas  for  the  series  of  asymptotic  forms  of  the  Beseel.  functions* 
The  difference  series  may  then  be  more  quickly  summed* 


The  formulas  for  2 and  the  various  aeries  involved  are  first  tabulated  and  put 
into  dimensionless  form*  The  asymptotic  forms  of  the  Bessel  functions  in  Coo  are 
introduced,  giving  a Fourier  series  as  an  approximation  to  J?.  A numerical  part  of  this 
series  is  sunmed  by  use  pf  the  'f''  function,  its  derivatives  and  asymptotic  form,  and 
ths  result  is  checked  by  a second  method,  using  the  Euier-Maclaurin  foraula*  The  Fourier 
series  is  summed  by  starting  with  simpler  series  with  easily  obtained  sums,  and  then 
building  up  the  desired  series  by  successive  integration,  leading  to  rapidly  convergent 
power  series*  Two  forms  of  the  power  series  are  given  with  different  ranges  of  convergence* 
Numerical  results  obtained  from  these  formulas  are  discussed,  tabulated  (in  Appendix  IV) 
and  plotted  for  a range  of  C from  1 to  10  and  J£'  /?'  from  0 to  / 6(/+±)  „ The 
values  are  intended  as  a survey  of  the  physical  behavior  and  as  a teat  of  the  formulas, 
hence  are  not  pushed  to  high  accuracy,  but  are  based  mainly  on  the  first  term  of  the 
variational  expression  for  2 depending  just  on  C&* , The  basic  terms  of  the  Bessel 
function  series  are  given,  from  which  the  first  order  and  all  higher  order  calculations 
wuld  start*  Sufficient  data  are  given  so  that  2 may  be  simply  calculated  to  first 
approximation  in  the  general  case  of  change  of  cross  section  and  material  (four  parameters)* 
(^/z.)  and  tfZ  (for  matched  termination)  are  themselves  tabulated  for  Just  the  changs 
of  cross  section  (two  parameters)  over  the  parameter  ranges  given  above*  The  second 
approximation  is  evaluated  for  various  values  of  C and  , and  makes  at  most 

a ZO%  correction  to  ( 2/z.)  (Z  % fo  &) 


Procedures  for  evalu  tion  of  the  C&r*  for  any  values  of  & and  /n  are  now 
considered,  since  these  must  bo  known  to  Owummi  who  higher  order  terms  in  2?  * Again 

the  asymptotic  (Fourier)  series  is  obtained  and  summed  by  procedures  similar  to  those  used 
on  Coo  c The  general  sum  for  Col  is  expressed  in  terms  of  the  ¥ function,  trigo- 
nometric functions  and  a definite  integral*  Two  procedures  for  calculating  the  latter 
are  given  in  Appendix  II,  namely  series  expansion  and  numerical  integration  by  Filon's 
formula*  Special  consideration  must  be  given  to  the  ofic\y?ence  of  an  accidental  singu- 
larity in  one  of  ths  terms  of  the  series  due  tc  the  use  of  the  asymptotic  forms;  it  is 
proved,  however,  that  the  exact  series  is  always  finite*  The  series  of  the  coefficients 
Ctt  is  sunmed  by  differentiating  the  series  for  Col  * The  quantities  Cjm  are  easily 
obtained  from  Cot  and  Com  by  a simple  formula*  It  is  noted  that  the  Hahn  functions 
are  easily  sraused  by  the  methods  of  this  section  and  formulas  for  them  are  given..  Finally 
the  additional  factor  in  the  sums  which  depends  on  frequency  is  discussed*  This  factor 
ia  easily  taken  account  of  numerically  when  the  wavelength  is  long  (compared  to  rod  radii), 
but  is  more  important  as  we  go  further  from  the  static  limit*  Summation  procedure  for 
series  containing  the  next  tern  in  the  expansion  of  this  faefor  is  indicated* 


6b)  Recapitulation  of  Formulas 


We  have 
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where  the  roots  with  positive  real  farts  are  taken  in  (6.7),  so  that  the  wavs  £>*  ^ * 

attenuates  in  the  prorogation  direction,  namely  toward 


Thus  the  first  approximation  to  ^ is 

" if  ,‘r*-£Dt.  ^ /D3  J -/  W 


(6S) 


2 ' a ***  f.a &j/( £„)[,- 

We  proceed  now  to  the  numerical  evaluation  of  (6.S),  which  is  a series  decreasing 
as  l/n  , hence  some  20  terms  could  be  required  for  4 place  accuracy. 

(6c)  Introduction  of  Asymptotic  Forms 

The  convergence  of  (6.8)can  be  improved,  however,  by  subtracting  the  asymptotic 
form  of  the  general  term  and  stunning  analytically  the  series  of  asymptotic  terms.  We 
have  the  asymptotic  expansions 

S -i'«  * <* 
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Taking  the  leading  terms  in  (6.9)  we  have 
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Hence 
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(6d)  Evaluation  of  Numerical  Series  by  Use  of  the  sP  Function. 
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(See  -Jahnke  and  Ende 
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(7)  E.  Jahnke  and  F.  Erode  "Tables  of  Functions,  3rd  Edition,  1933 
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. (6f)  Summation  of  Fourier  Series  of  Asymptotic  Terms,  First  Form  for_.  0 

» j-//?  lir/rff.  '/*  ) 

The  evaluation  of  the  second  term  in(6.12),  namely  1 


/9=/ 


requires  more  elaborate  analysis.  We  note  that  if 

f S/n(n+'A)} 

£i  (fl+'A?  J 


(r>+  'A)' 


then 


/ £ cos(fi+A)£ 

**'  (n*'Mz 

V<»'  -f  s’inlnr'A)^rU) 


fiH 


/"()) 5 ~2  cas(n+‘/+)j  - -cos  ) + Ssn  k(2/  Sfn n*) 


(£.20) 

(£2  0 
(£22) 

(€.23) 


The  series  in  (6,23)  are  not  convergent  in  the  simple  definition  of  convergence  but 
rather  oscillate.  However,  they  may  be  summed  by  a slight  extension  of  the  usual 
definition,  and  then  integrated,  yielding  in  (6,22),  (6.21)  and  (6.20)  series  con- 
vergent in  the  ordinary  sense. 


Thus 


jV7^*:  -fp  = T + 
£ 0 /-*  * Z 


~ cot  3L 

2 2. 


or*  formally,  we  may  write  (using  the  symbol  ( ~) 

OO  0*3 

'Z.cosrij.(*)f  ) 2 sinnj-(*)-zcof 2 

“ /Ji.S  ^ ■’“* 


(6.2+) 

to  mean  ft  is  formally  equal  ton) 

(£.2  S’) 


n-o 


n-o 


Mote  that  (6.24)is  convergent  if  a has  a small  positive  imaginary  part,  although  this 
is  not  true  for  (6.25).  The  meaning  of  (6.25)  will  not  be  made  precise  here,  but  we 
shall  merely  use  the  formal  equalities  to  suggest  new  equalities  obtained  by  operations 
such  as  integration.  The  new  equalities  oan  then  be  checked  and  established  independently. 
Regardless  of  the  way  (6.25)  is  given  a lire  else  meaning,  che  equality  clearly  cannot  hold 
4t  3--0  or  2 =2?r  and  may  be  expected  to  hold  only  over  . Then  putting  (6.25) 

in  (6.23)  will  give,  in  the  sane  formal  sense,  (note  the  change  of  sign  and  summation  limit), 

£ cos(fJ+/+)2(=)£cos l - S>nk(±cc£ k)  s ( !^co. t) ; 0<}<Zir  . (6,26) 


n*<> 


Putting  t - 2/4. 
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As  a check  we  take  ~K  - cos.(+n+/j^ c)n fz/i  , ond  now  w&  must  hcn/e 


A..{lz 

Z /f-o  f +0+0 


TTZ  j_  ^ £„  i 7Tl 

3 2 r + £b  (2n+Z)f  ( /*  / 


5 f["£ Ti  +£'• 126,7  ” /! " 7** 


which  checks  (6.3=) . Thus  we  have  reason  to  believe  (6.35)  with  (6.36)  for  A 
is  orreet  equality. 

nally  integrating  (6.35)  from  0 lo  x,/not-«  (6.35)  now  converges  at  gives 


f {&J*ZLtil3L  -/o  7P8  3 Z -J-J  _£>»  * 

/'  •«  _ . * i j?  / <$  a ^ /i  ..?)y 
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(Znt3)/ 


To  checks 7^,  toAe  Z*  *5$,  or,  J pur  S//)ten+/J7&  - f-Jn!^/2  , 
so  that  the  sum  becomes 
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z nTo  (+nr*)*  2 !/»•« 
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(6.37J 


S/2. 


= -7020 


(6.  3S) 


to  be  compared  with  loO?/*83x  .78540  - Q19635x  .61685  - 1/4  x .08344 
« .70219  directly  from  (6U37) . The  coefficients  in  (6.3?) 

En/(  2/7/-J  )/  are  tabulated  below  ( and  in  Appendix  IV) . 

n ^>4znr3J/ 

0 .166  666  667 

1 .008  333  333 

2 oOOO  992  064 

3 oOOO  168  100 

4 o000  034  697 

5 ,000  008  113 

6 ,000  002  067 

? o000  000  561 
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(6g)  Second  form  of  Sun^jdon  of  Fourier  Series  for  Near  / 

The  aeries  (6.37)  for  § (the  asymptotic  series  to  S)  does  not  converge  well  for 
77/£/"/X  or  /" <7^  z /.S~70&  (thus  the  eighth  term  is  /r7,-7'('%)/T  = - 

To  cover  the  range  / <-  /'i  /r/2  another  aeries  is  useful  which  may  bs  derived 

similarly,  starting  from  (6.27), as  follows; 

I m ; £5* . X *i  l v*"'/  y - f-*  • ^ ■** 

/)to  ~ s*z/ 

Integrating  from£«^(  ^ or  frcmy- ^ /» y as  in  (6,30),  and  using  (6,32) 
r"  JV/7  fa/??/ J 2.  /p  /"  7T  / i 7 ) t y i/*!  & ^ 

£ -x^7“  - itUyt  *• 

¥ allil  £k  u*"\ 
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(6.40/ 


Note  that  /v  ^ - /ny  - /nZ-r^Zz.  — ; , hence  the  two  forms  of  (6,10)  give 

j£//i ( fi ~/j + /o 2 - />?</ - ”-/ny -%'-■]  (6'+/J 


Thus  the  numerical  series  in  (6.11)  or  (6.40)  is  evaluated  and  (6.40)  becomes,  after 
cancellations, 
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Integrate  (6.42)  from  2 to  or  from  v’  fv  0 


at  y»0  is  integrable,  to  give 
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noting  that-  the  singularity 
2/7  H 
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Integrating  (6.43)  over  the  same  ranges,  uoing  (6.1?)  gives 
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and  compare  with  (6.38),  giving  in  both  cases. 


To  check  (6.45),  Lake  ‘Z-y  * 77/ 4. , 

O,  -702  /9. 

The  coefficients  in  the  series  in  (6,45  ) a re  as  follows : (also  given  An  Appendix  IV) 

Bn 


*- 

3 

4 

5 
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(in  *-*-)! 

,006  944  444 

162  037 
006  102 
000  292 
000  016 
000  001 


=000 

„000 

„000 

=000 

=000 


These  are  substantially 
smaller  than  the  coefficients  in  (6.37), 
a result  probably  due  to  explicit 
removal  of  the  logarithmic  singularity 
at  r ^ or  y*0  in  (6.44X 


Both  series  (6,37)  and  (6.46)  converge  over  the  entire  physical  range  of  <jT  , but 
(6.37)  diverges  at/fc*7^  or  /"»/,  while  {6,46)  diverges  only  when  y reaches  77“  ; the 
largest  value  of  physical  interest  is  y - at  /^=  00  . 

Thus  (6,46)  is  more  rapidly  convergent  for  /!</" i 3/(oi y $(”7z))  while  (6,37)  is 

bet te r f or  J 4 f £ (V‘ ><?,). 


(£h)  Evaluation  of  i?  and  (K,  General  Formulas  and  Physical  Parameters 

We  can  now  proceed  with  the  calculation  of  Z , or  better,  of  ^/z»  , 

From  (6.8),  and  (5.15)  - ^ , 

7 /r  f ^ (*n/r) 

» -*KCoo  A,  17^* 

r»  ^ ; J,(Ay>o 


•//— 


Thus  ^2^  depends  on  four  parameters,  (%■■),  r,  jt/?,  jcr' 

%,  = -****'('&■)>*■>  * if, £? ‘I Snl/^A)^  ‘ 

Now  is  asymptotic  t-o  £ , given  by  =5“  =/^/  J C&*))  - ^ 
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From  (6-453#  however, 

/- 2 /)  2 j * 2/9^2 

>■*/-{-,. a3i°*sO-  ™ ?-)*■  -.  J — _/  (fcjv,’ 


wh«vo  yf  ' 7^2/"  s ^ ~ ^ • 

v 

Thus  jj£  depends  on  just  the  two  perimeters  <T~  ond  Jt  ‘/x* ' , while  •S'^v  -5* 

depend  only  on  /*  . This  makes  it  possible  to  tabulate  x*  as  a function  of  </''  and 

it'/?’  and  permits  evaluation  of  2/2*  and  the  reflection  coefficient  for  all  changes 
of  cross  section  and  of  material  by  a simple  calculation  (illustrated  in  subsection  6i) 
which  brings  in  the  two  additional  parameters  K/K*  and  SC  R-  However,  if  there  is  no 
change  of  material  but  only  a change  of  cross  section,  then  >d  ■ t X«K’,  and  z/20  depends 
only  on  S'  and  Jt '/?  ' (cr  dLK  ) and  the  reflection  coefficient  may  be  tabulated 

directly  as  a function  of  these  two  parameters  (Appendix  IV)  « 

These  considerations  make  use  thus  far  of  just  the  first  approximation  (6-46) 
to  */cm  in  the  variational  form-  The  effects  of  higher  order  terms  will  be  considered  later, 
but  in  general  will  amount  to  small  corrections* 

The  voltage  (or  temperature)  reflection  coefficient  for  semi-infinite  termination 
of  the  rod  /?'  is  then  given  by  (4.19),  (4.20)  or  (4.23),  and  takes  the  fonns^for  the 
special  case  of  a single  series  element  2 at  the  junction, 

/+&  Z . * - * J£/?L  (d>.o-£j 

JTZT  & 2o  ~ & /T  ■ 


We  note  the  second  term  is  real  and  for  a change  of  cross  section  without  change  of 
material  is  Just  ( ' /tr*- ) 


(6:L)  Numerical  Results 


Sane  calculations  have  been  made, based  on  the  formulas  of  this  Sect ion, in  order  to 
survey  the  behavior  of  the  series  impedance  representing  the  junction  over  a reasonable 
range  of  the  physical  parameters.  In  addition  the  calculations  serve  to  test  the  usefulness 
of  the  formulas  for  numerical  purposes.  In  view  of  the  large  amount  of  numerical  work 
required  for  a complete  and  accurate  tabulation  of  the  impedance  as  a function  of  four 
parameters  ir.  the  general  case,  and  since  this  is  intended  as  a survey  of  the  general 
features  of  the  problem,  the  work  has  not  been  pushed  to  high  accuracy-  nor  has  it  been 
careful  y chocked.  However,  enough  has  been  done  both  to  give  these  general  features  and 
to  serve  as  a starting  point  for  an  accurate  result  for  any  particular  set  of  parameters. 

Accordingly,  most  attention  has  been  given  to  calculations  using  just  the  first  term 
of  the  variational  expression  for  ? , as  in  (6.46),  but  at  a few  points  the  magnitude 

of  the  second  term,  aa  in  (6.1),  has  been  found, to  give  sene  idea  of  the  accuracy  of  the 
first  approximation.  In  Appendix  IV  a tabulation  is  given  of  the  terms  S#  (fj  defined  in 
(6^48),  over  the  range  of  d'  from  1 to  10,  and  also  of  the  terms  of  the  asymptotic  series 
Cc  of  (6.49),  to  illustrate  the  summation  procedure.  The  Sn(f)  are  the  basic  quantities 
entering  all  the  numerical  aeries  for  thy  coefficients  Cam  and  the  starting 
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point  for  any  more  elaborate  or  precise  calculation.  From  them  the  suma^f-^ 

(or  ^c)( CA'ft’) ) are  obtained  for  the  same 

range  of  /"  values,  and  for  a range  <s  £ from  0 to  1.6  ( UrJ)  corresponding 

to  a wave  length  comparable  to  the  rod  diameter.  j^Since  fSS-D 

the  maximum  value  of  corresponds  to  a maxisaaa  frequency  of  , which  at 

low  temDeratures,  for  metals,  and  for  reasonably  sized  rods,  is  at  least  a fsw  thousand 
cycles  per  second.  This  is  about  as  high  a frequency  as  the  thermometers  in  use  can 
respond  to  effectively,  although  further  progress  here  might  extend  the  ran.-e  of  interest 
of 

Fran  (/(Jt'ft'j  , by  (6.46),  ^ r? / ' ) can  be  easily  found 

for  any  values  of  the  four  variables.  If,  however,  we  take  the  special  case  in  which 
no  change  of  material  occurs,  then  &ft  • JCft'/f  and  ( */A'J  * 4,  and  (*?£,) 

becomes  a function  of  just  two  variables.  (*¥2*)  is  tabulated  in  the  next  table 
for  this  case,  over  the  same  range  of  /"  and  X ‘ft ' as  used  for  J,0'c,j 

and  also  the  corresponding  reflection  coefficient,  & , for  a wave  incident  on  the 

junction  from  the  smaller  rod,  when  the  larger  rod  has  a matched  termination.  The 
results  aj*e  plotted  and  show  that  ( ^/z.)  varies  from  o of  S*/  , through  a mar! mum 

at  about  i'-Z.S'  , to  smaller  values  as  /"  increases  to  10.  2 remains  always  in  the 

first  quadrant,  as  required  by  the  theorem  of  Section  4,  and  in  fact  is  nearly  real.^ 

(2/2»)  is  quite  frequency  dependent,  rising  from  O at  aero  frequency.  However  ,,>J 
is  much  less  frequency  sensitive,  as  shewn  in  the  plot,  and  varies  by  less  than  3 1/3& 
oyer  to  1.6  (/r*  ) (in  the  real  part-  which  is  th£  dominantpart). 

The  values  of  rise  fr<m  o at  ! rapidly  to  nearly  *"*/  as  / increases  to  /o 

(the  large  value  of  / corresponds  to  a short  circuit,  i.e„  ai  end  held  at  fixed  7~ 

by  the  large  mass  of  the  secerid  rod),  and  also  varies  substantially  with  frequency  in  both 
magnitude  (which  decreases  at  higher  frequencies)  and  phase.  Thus  at  <Tm  2.5,  the  reflection 
at  ■ / C( l+i)  is  down  from  .7  to  .4*  and  varies  nearly  linearly 

with  W 

Finally  the  second  approximation  to  is  obtained  for  selected  values  of 

f and  Jt '/?'  . This  requires  tabulation  of  ( S',  Jt'R)  and 

as  defined  in  (6.54)  end  ( "Z/2*)  can  then  be  calculated  for  any  values  of  ths  four 
parameters  to  second  approximation.  The  appropriate  fo;rmula,  obtained  by  transformation 
of  (6.1),  using  (6,53).*  (6.54)  is 

(pfo.O^ 

I ■ -*  =r  / 

2 . A L J + i 'A  J 


In  the  came  special  case  a:;  above,  i.e.  no  change  of  material,  ( ¥2*)  ana  JC 

are  tabulated  explicitly  for  tha  selected  values  of  S'  and  X ‘ft ' . The  results 

show  the  second  term  in  ( */**)  is  at  most  about  10 % cf  the  first  tern,  being: 

emails r at  larger  values  of  f and  lower  frequencies.  In  /^/  the  corrections  are  at 
mist  only  about  2£« 


‘ v 


F/g  6-2  The  Comp/ex  Bessel  Function 
Sum  S(*oi( F,  H' F'J  Which  Enters  the  Ist 
Variational  Approximation  to  Z . 
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Fig.  6 - 5.  AAagn/tude  of 
Ref  feet  ion  Coefficient  in 
Same  Situation  as  F/y.  64 
as  a Function  of  Frequency 
for  Various  FT 
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(6j)  Higher  Order  Coefficients  in  the  Mode  Amplitude  Equations 

The  general  coefficient,  in  equations  (5»31),  , must  be  evaluated  to  obtain 

higher  mode  amplitudes  and  higher  order  approximations  to  the  variational  expression 
for  J?4(6.1).  The  general  formulas  from  (6 >2)  and  (6*3),  are 

c,m ' Xs uw  f-  a*  ](~^)[, 
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Note  that  if  the  new  factors  in  the  denominator  (of  Srf^'/n/  compai^d  to  5o  ) 
vanish  or  beeomB  small,  the  corresponding  term  in  is  still  not  singular  because 

the  mime  rater  -"’rishes  o?  5-  ~zAi2  » Thus  if 


then 


,/  )•  U(&  * (fyV=  J,  (p£)  +J,  ’(&*•  %)  = ^ ^ ^ ^ 


(l,a) 

and  the  term  in  5 is: 

£2  J>l(M  /_ 


£ /?£  Jo  "Ufa) 


W j *■(£„)  _££ /, _ CjQzL  ) ^nJ**^*)  (/~  ^jfp  ) 

r A Jo  x mP  " 

which  is  small  of  order  € „ The  cor  e 3 ponding  term  in  s will  not  be  small, 

however,  but  it  will  be  finite. 


(6.£Sj 
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(6k)  Summing  the  Asymptotic  Series  for  the  Higher  'Order  Coefficients,,  Com. 

Wo  give  now  a procedure  for  summing  l ha  »sf ‘pletic  serisc  for  the  higher  order 
coefficients.  The  formulas  are  somewhat  complicated,  and  would  be  most  useful  when 
th  •>  corrections  to  the  first  term  of  the  variational  expression  for  2 are  no  longer  small, 
hftncr  must  be  accurately  evaluated.  This  occurs  when  the  wavelength  becomes  comparable 
to  or  smaller  than  a rod  diameter  and  several  modes  are  excited. 

These  formulas  are  also  of  interest  for  two  other  reasons.  First,  they  provide  exact 
sum  formulas  for  the  coefficients  of  the  mode  amplitude  equations  in  the  case  of  the  change 
of  section  of  a cylindrical  rod  of  rectangular  cross  section,  where  an  exactly  analogous 
development  is  possible ( the  change  of  section  must  take  place  in  only  one  of  the  rectangu- 
lar dimensions).  It  should  be  remarked,  however,  that  more  powerful  approaches  are  possible 
in  this  case,  namely  the  solution  by  the  equivalent  static  method  and  conformal  mapping. 
Second,  the  formulas  provide  exact  sums  for  the  Hahn  functions,  which  vere  conveniently 
introduced  in  the  solution  of  a cavity  resonator  problem  by  mode  matching  techniques. 

These  functions  are  defined  later  and  the  appropriate  sum  formulae  are  given. 


For  definiteness  we  take  the  coefficient  C0> 
in  the  form 

Co,  _ _ C ft-s;n(n+'/+)j) 


t£  R:Joi/J 
/r 


- /<■>, 0 . ZJT 


and  write  the  asymptotic  series 

(S6  oj 


A proce  *e  for  summation  of  thi3  series 

>(*,0  . r /'r  L f 

3 n--/ 

where  XT  - (*'^77j  is  approximately  /4-) 


is  the  following.  We  have 

Sin  (s?r  VW>  ' . » 

(n+'/<y/7t'U  +/£)(n*~/4  ~£) 

(we  shall  use  the  notation 


and  f £ - 2.^, . 


Now 


“ 2jC  *1  ' m'A-f  W 


fm 'U)(r>t '/4+/: tfnt'U  -/>'  2/C 


(6,  6 2) 


hence,  using  (6.1A) 

/ £,[-  m r,c)  - m -t)t2fr/')] . 


GO 


(n+V+Xn  t Z*  +/:)(n+  # ~£) 


(6.63J 


Again  using  (6.61)  the  other  term  in  (6.61)  is 


V*  5 /V»  (0+ 

Cnt -A'/zI^/a T/r)(nt'/+  ycj 
/?•-/  ' “ ' ' 
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The  last  tern  in  (6.64)  has  already  been  evaluated  and  the  two  expansions  for  it 
in  (6.34)  and  (6.42;  are 

G(i.» E f ^ ”"J 

•-  [f  - &7  <‘-‘V 


where  the  new  function  of  two  variables  G(a.z)  has  been  introduced. 
The  other  two  terras  in  (6.64)  are  of  the  form 
% . cos Ci  -a)2?  -ok  ~ 


r>-/ 


f)+  4 


/?=/ 


/-  / 


= rosfi  -a)£  G>'oJj)  - T/r/ ^ 'crt'a'j) 


(6.  6 6) 


vrb«»**  •nether  function  or  erro  variables,  F(a,z),  is  introduced . 

O Z C.*y+  , /u  - C - fas 


Put 


and  we  have 


^ Sir.'ft+'UJZ  . <T  Dr  G f~  Ajjj  - S'/O/:  J f %J)  - )j 

/L  — — - , ’ T~  s>  + 'S+  • <r  . *"  ' 

/,=/  vt’U+F  nx/  ’* 


The  sums  F and  G may  be  evaluated  as  follows.  From  (6.25) 

jf  cos  C,-?+*JZ  (-)  fCCos  -SJnoj- 

f s/* f&ojj.  -rstnai) 

/j*o 


(e.eaj 


and  upon  integration  from  7T to  e (the  series  (6.68)  do  not” converge”  at  a»0 


or  2rr  > 


<^,*22  ?z  aj-  / (e 

9 nio  ' r " O \ 

F('a,2)+CC6.aJ'z  £ f ‘ coson2  .if  < coso^co^  efr  ^ 
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The  first  terms  in  (6.69)  are  given  by 


( — v 
flta 


, n o* 


°*  / ®£ 
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m>/ 


sr>'/ 


2.mro-/ 


- -t  ?(¥h  n°-fj 


(6.  ?o) 


The  integrals  in(6.69)  are  put  in  a form  for  numeric&levalvtatiou  by  two  different 
methods  in  Appendix  II.  One  procedure  expands  cotJ’/z  in  pavers  of  ^ , and  inte- 

grates term  by  term,  obtaining  a rapidly  convergent  series  over  the  range  of  interest 
of  ) t The  other  method  sets  up  & numerical  integration  scheme  specially  adapted 
for  integrands  with  trigonometric  factors  s.-n&j.  and  eosoj~  . The  method  (due 
to  Filon)  is  as  simple  and  accurate  as  Simpson's  nile  applied  to  the  integrand  without 
the  oscillating  factor. 

In  both  cases  it  is  advisable  to  subtract  trie  singular  part  of  cot  J/z  namely 

and  evaluate  this  term  separately  in  terns  of  cosine  and  sine  integrals. 

This  gives 

L fe  i0^cot<£  Jj,  = - c*  (Jf})+  i ( S<( ai}-Si(aTT))  rje  *(j  **  2 cc7t  \)di  ( 6.  ~7/) 

where  the  integrand  ^ tot-*/z  is  a finite,  smooth,  monctonic  function 

of  7 ovsr  o to  77  ranging  from  O at  ^ = to  .32  at  j - nr  0 It  is  shewn  by 
jor./  that  only  values  of  3-  between  v and  7~r  need  be  considered  in  (6.71) 


The  srauns  G(-bj) 
substituted  for*  a. 


Then 


6(-A})‘ 


and  ?(-&,£)  are  given  by  (6.69)  with  — b 

_ sinbrrf.)  ?(—-)] -■*—)  3 (6- 


1 “‘4  «*»  1 1 the 


functions  be  removed  by  us©  of  the  relation 


\p(~  x)=  '£'(%’-/)+  Trcotnx.  - zrz?  /'Tcot/rx. 


Hence 


-it (=&■)+*  -itr-rH 


Thus,  from  (6.69),  (6.70) 


and  using  (6.74) 

✓ , sin  62. 

= - G,%}  )rrr  - 

Similarly 


F(-^)  - r:°3 


) _ . fi  jV  - TTCutjt’77 


((  73} 


&W*  s>n*n(.i  W))  ~f  U.  7rJ 


(6.  74) 
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Collecting  terms  from  (6064)  (6*65) * (6067)*  (6.75)  to  (6=78) 

f - *loJjl±!F±LL * -±  £/-'  j)  y -A  c*s/=j (g fa; )-&*>> 

{n+'U Xn+ '/+  t-/c)(n  + '/+  -/c)  /:*  ( 6 / 

S//tCj,^fotj)-rt4j)  +7TC&f£>7r  - 


(6  79/ 


whereas 


? ^ ■; , 6 /"  J>(Qjr&hJ--l*-77Cert>rr  -l  ?f£)7  (g.&oj 

f-,  (ntWAtytr/cXor*-/:)  Z/d[  r r ~ +'J 


r>-/ 


and. 


7 f r».£o)  - (6.  79jJ. 

£tt*L  • ' 


(6.  eo 


(6.79)  to  (6.81)  with  (6.77),  (6.79)  and  the  procedures  of  Appendix  II  permit 
evaluation  of  the  asymptotic  s«rie3  for  Co/  , or  more  generally  for  Com  , Just 
as  (o0i/9)  provides  for  Coo  * 

* 

(6£)  Difficulties  for  b Nearly  Integral 

We  note  one  further  complication.  If  by  accident  b should  approximate  an  integer, 
cot/>/7-  diverges.  Now  the  two  terms  in  (6.79)  and  (6.60)  containing  cot bvr  do  not  cancel  if 
sin  os  is  not  unity,  and  the  asymptotic  series  is  .then  divergent.  However,  it  was  proved 
in  (6.59)  that  the  terms  of  the  original  series  S'0*1',  are  always  finite,  so  that  the 
error  comes  in  through  the  use  of  the  asymptotic  forms.  Accordingly  ths  asymptotic  series 
should  be  corrected  around  this  sensitive  term  by  removing  these  terms  and  putting  in 

•*  . <*_  fcnr  odrilni?  fcr»  the  SITU 


correct  ones,  i.  e.  by  adding  to 

j, > 


y l J (&•'*')  — j£  OjJ  J 

C,i  (aV.  Zn'rm'AIfnr'brtJInr  * -*)  J 


(6  82) 


A3  many  terms  of  (6.82)  should  be  added,  for  n above  and  below  b,  as  contribute  in  the 
decimal  place  of  interest.  Note  that  if  b ■ c-l/4  the  term  n**£  in 

is  approximately  <^#3  which  on  subtraction  from  the  sum  formula 

cancels  the  sum  of  the  two  cot  terms  in  (6.79),  (6.80)  Thus  any  singularity  is  removed. 

f Oti J 

If  we  make  the  further  approximation  in  S of  using  the  asymptotio  form 

of  c,  namely  c - » then  yCCKi^2  =<3  s/s?Cj-  and  (6.79)  to 

(6.81)  give 

'-«»)_  r « /-SinC**'/*)* 

^ ~ zn*  /Ft  (n+ ^ " ‘S+) 

,jCrri[-f(S/-*(i)t2Wi)*  (6.83) 
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where/T-  3 = 27^. 


In  (6, ,33) the  possible 


singularities  at  b • an  integer  have  cancelled  out,  but  at  the  cost  of  a further 
approximation  which  may  be  a non-negligible  fraction  of  each  term  of  the  series 
(actually  a few  percent  since  c*1.220/" ✓r.  j?  -/zs’cr'  ).  The  fraction 

decreases  with  increasing  n,  however,  once  the  critical  term  (/?■="*  ) has  been 

passed,  so  that  for  small  f (less  than  three  or  four  say),  (6.83)  might  je  quite 
usable  as  an  asymptotic  series  for_S  or  . A numerical  example 

applying  (6.83)  to  calculation  of  is  given  in  ; ppendix  III. 


(6m)  Evaluation  of  „ Cjtt 


Application  of  this  above  procedvire  to  evaluation  of  more  general  Cjln*  will 
prove  very  simple  when  1 / m,  but  more  difficult  for  Cub.  Consider  the  latter  first  > 
and  in  particular  the  prototype  case  C0  . Fro&v(6.57)  the  asymptotic  series 

for  C7/  is 

c„  -<>,>)  r ~ (n+  '/+&)( /)t  fAJ 

4-Z.  ~ ~ £:,(/*■> ’trAffr*'/*-*)1 

K‘  w S?  i i in  f/>*  •/*)!) 


(6.  S*J 
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err*/: 
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where 
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fl-i)  * / S3*  +s7&rr-  * - * ‘W  * s£dr 


n*t 


(6.88) 


and 


cos  (01-0)3 


<5,  WZ  ) 6 


0*1 


sm  (01-0)3. 

(nro)*- 

Integrating  (6.77)  from  0 to  z yields 

r , 7. /-!-■>■  ■ r/o-t  ) L 7 • sinaj. 

GJO^)  * ) cosorrj  - z r / ~ *(  irJ'ZoJ  ~Z  “ * 


(S.  69) 
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The  first  integral  in  (6.90)  has  been  given  in  (II  .10) 
we  obtain 


By  a similar  procedure 


yfcosoj-J.  cofts  d,  r y-£  Q-J±>(s*oj.  cosoj  S^Jvjj)  (d  Wj 

cj  & /?*,  cr” 


A numerical  integration  procedure  can  also  be  set  up. 

Integrating  (6.75)  from  0 to  z gives 

F2  (<zj) ffaj  'F  ^a7T[-{  ‘S' (^r)  + T T”7  ~£]~  ~~2~£f  i 

/■  ftiooj.  cotjc.  e/j.  - j/ A cob%  *0  fcfZJ 

Z Jr  * o 

where  the  first  integral  is  given  by  (II  oil)  and  the  second  integral  by 
y oa  r> & 

= <Z^i££  -2  (e,93) 

CJ6  a />-/  a 


"he  formulas  for  F(  ~ At ) A -F  are  then  obtained  by  integrating 

(6.78),  (6.76)  respectively.  This  gives  immediately 


6;  (’Aj)  = Gz  (6,})+  s_F?bJ  - 77 £ coStrr 

bi 


and,  using  (6.83) 


^ ^f^j)  - - —p^-  t jfn^Z/r 


(<$. 


(d.  s>sy 


Finally  putting  (6.88)  (6.94)  and  (6.95)  into  (6.86)  we  have 

gO,,)  /-  f’M- £ r (*.;)- 


(6,96) 


Note  that  in  (6.96)  for  the  special  case  /£J  - ~ , the  singularity  when  b 
is  an  integer  cancels  out  as  in  (6.83).  However,  in  the  crse  of  general /£ 2.  , we 
have  the  singularity 


^-[F^r^'S'r,SjA  cos;;  frrj.  ce+A*)J 


&77J/C 
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v'y*sVl 
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<r  (/-*">( £t '/+)j>)  / 

ir  the  termr^j  77  772  in 

ifr  (j£+/+ 

is  subtracted  from  the  atm,  -his  subtracts  the  singular  t 


(6.85),  for  which 
r r/s^j-  z co^j  7 

Vi  IliU  — > J ■“ ■ — — __  j 

8rfcL  €k  £ J 


from  the  same  term  in  (6.97),  and  removes  the  singularity.  As  in  (6.82),  the  correct 
Bessel  function  terms  can  be  used  in  place  of  the  asymptotic  form,  for  the  terms  near  £. 

(6n)  Evaluation  of  Pin  /t? 


The  calculation  of  Zem  can  be  made  to  depend  in  a direct, 
and  <Z0ni  , so  that  it  is  unnecessary  t,o  sum  the  asymptotic  series 


transforming  (6.53)  by  partial  fractions, 

r tA  J?U0(/?e)UlUm)  1 ^T~~Y  fa 

^ /r  Aj  JUM  /?„(*/ -as)  L& 


simple  way  on  C o£ 
. We  have,  on 


A A 7 


- A ~~  z fj%t  A (&*•)  Col  J-  (At)  Com  f £ p 

L */  j 


(s.ssj 


Thus  the  calculation  of  the  two  dimensional  array  of  coefficients  C in  the 

linear  equations  (5.31),  is  reduced  essentially  to  the  calculation  of  two  linear 
sequences  of  coefficients  C5/  , and  C n 

(60)  The  Hahn  Functions 


We  note  tliat  the  procedure  introduced  in  the  last  sections  to  sum  the  Fourier 
3ories  asymptotic  to  C/^  can  also  be  applied  to  sum  t e Hahn  functions  introduced 
scasr  yearn  ago  by  W.  Hahn  in  the  solui-iQn  of  a cavity  resonator  problem  by  mode  ex- 
pansions and  field  matching  methods, ' ' 

Hahn's  functions  are  defined  as 


S,(°)  - 2 


Sin 1 rrma 


alsn  S/nrrma 

; u (°)  * / Ti/jjZEF  z7  ^ ***9"  (c 


PL  / ' ( >*  v 

and  in  the  reference  a sum  formula  for  Ca)  is  given  but  not  one  for  Up  f a)  ■ 
B(y  analysis  precisely  like  (6.61)  - (6.64)  we  have 


Sp  (*)  * / L 0 . - casmj.)f. 

//? ' / t 


r>itcr 


/ __  z 

/77-  o'  /r> 


zi  tei2/ro 

/ <r  --  r/a 


( 6.  /oo) 


and  using  (6.14),  (6.73).  we  have 

j- 1 (sk  -k)  - tU  tlrl+t-rrcctmri.  - pfotJ 
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4 

i 


while  analysis  like  (6.66),  (6.67),  with  (6.76)  and  (6.78)  gives 

/ 1 ' *rht  - -kj*  (zw*  - srcorvTTj 

mu 

y-  Sir.  ^ (z  3.)-7Tt  £^Lf^  ) . fa  /0*) 


(6.101)  and  (6.102)  in  (6.100)  give  the  sum  for  general  values  of  cr  and  . However, 
in  the  original  definition  of  (oj  p is  taken  to  be  an  integer,  which  leads  to  a 
simplification.  Thus 

O'?.  = 2 77/?,  COS  tr£  - /,  s/f)  <rjr~0  (£•  703) 


«|na 


Sp(a)-.-i^rcr)-jf^p  -£/mZ</-cosj)-  - 
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(&/04J 


- ^(p/j)  ■ i /-( %)Zno)  - £/n[z( /-coil 7To)J - 

3-  *"  ^ ^ //?  0-/7  integer) 

where  1 7To)  is  evaluated  by  (6.77),  and  (11,10  ),  and  $f°J  5 "<Csr  --,f77c  /6  has 

been  used 


We  note  the  special  case  a»l/2  for  which 


. c? 

"T  y 


Now 


F (2 A *)*[“{  /V>  -0+  / p (p -'/>)]  - jy, 


hence 


5*  = - t -M  Z-  ; />  or>  /sif  ter’ 

S'  U £ 

a result  which  is  readily  checked  directly  from  (6.99). 

Similarly,  we  have  , using  partial  fractions  just  as  in  (6.85) 


(6.  SOS? 
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which  becomes,  using  (6.87),  and  a formula  like  (6067) 


a f'f-o-)  -pftr)  ~/5  4v jj)  /*  S/n  <rj,(6z  (-<r,j)  + £. 


and  using  (6.88),  (6.95) » and  (6o94)  to  remove  the  negative  arguments,  gives 

'J  M -■  f[-2  f ^ -J,  rJgL.  /•  cos^( ?/f  (a#  rrrj.  e <2££Z  - 


ys,n<r£  fz  6zfof)<-  •r-^>  - fT£  cof-cr)J 


(6.107) 


(6.107)  holds  for  arbitrary  O’  and  2y  as  in  the  case  of  (6.102)  , but  for  the  special 
values  (6.103),  it  simplifies  to 


/Jp(d  = ?'(%  >)  ; 


/p  an  integer  (6.108) 


provided  thati?'s/%  is  not  an  integer.  VJhen/^is  an  integer,  the  term3  containing 


/-  co**-3. 

Gjn i c rrr 


and  3Jnrri.ca+<j-77  make  additional  finite  contributions,  to  give 


0,M' 


% integral,  (6.10?) 


f \ **  — yy 

(6p)  The  Factors  f / ~~  //£n)  ] 

The  calculation  of  ? requires  summation  of  the  various  series-*^  ' ^ in  which  y 
the  terms  S,i  are  multiplied  by  the  factor  involving  the  wavelength t[ t ~ J 4m. 

For  long  wavelengths  (compared  to  the  diameter)  M ’/? ,<c/  <cj8„  , and  the  factor 

quickly  approaches  its  asymptotic  value  of  1.  But  for  short  wavelengths*  this 
approach  may  be  too  slow  for  practicable  summation,  and  it  will  be  useful  to  subtract 

the  next  term  in  the  approach  tc  1,  namely  1 2‘  , and  to  sum  the  corres- 

ponding series.  These  series,  which  arise  from  3 and  o ' rn'  multiplied  by  this  addi- 
tional factor,  can  readily  be  summed  sines  we  merely  require  summation  of  such  series 
as 

y (i-sin(n*»U)2)  ^ (t-sifUn+y+te)  % (/-  S inf  At  '/+)’) 

Jr,  ' J-fin-kylintW- f ] ' J)  <n, 
flee  (GJOj  (QJ2),  (6. 60,  '684)) 
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The  first  aeries  is  found  by  further  integration  of  the  expansions  for 


r> 


4tlL_ 

7/  r 


ZL 


"he  second  requires  use  of  the  partial  fraction  expansion 

/ _Z_  / / __  ' 2 


1 ) 
/^V 


farfrWfi+'AJ1-**)  v+'U-*  Ot> 

and  all  the  terms  in  this  summation  have  beer,  done  previously.  The  third  requires 

/ „ - J_  f_2__  _ __/ L ' / ' — - L — ) 
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all  terms  of  which  have  been  done  previously.  Finally,  the  sum  of  the  series  with  the 
new  factor,  '/t  (*  - 


, multiplied  into  the  terns  of  3 
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immediately  from  the  sum 3 of  the  corresponding  series  obtained  from  3 ‘ ' and  3 {e* 
by  multiplication  by  the  new  factor,  just,  as  S ^ follows  from  3 <c’  ~ni  and 
(as  in  (6w98)j  . 


Then  we  have 


i s'.y^n  (*g'j  y 


where  the  last  series  is  eval- 


uated as  above.  The  difference  saviec  (between  left,  and  right)  will  now  Converge  mure 
rapidly  and  permit  easy  numerical  summation  for  larger  values  of  ( £ '/?'  ) - 

(6q)  Supplementary  Note  on  Evaluation  of  F(cf,j») , G(g~,>)  and  the  Hahn  Functions 
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After  completion  of  the  manuscript  the  author’s  attention  was  directed  by 
J„  R,  Vhinnery  to  the  work  of  Strachey  and  ’Wallis  on  the  Hahn  functions  A , 

By  methods  similar  to  those  above,  using  the  f function  in  summation  of  series, 
these  authors  evaluate  the  Hahn  functions.  Closed  expressions  are  obtained  whsn  ' O ' 
is  rational,  a situation  not.  considered  here.  For  general  & and  /£>  similar  but  not 
identical  formulas  to  those  above  are  obtained,  which  also  involve  a definite  integrals 
Tnese  may  be  shown  eouivalent  to  these  given  here  by  using  different  formulas  for 
F(<r,p  , 6 (<r,j)  , F2(<r,^)  r G 2(<r,^)  in  (6.10?)  , (6.104), 


Thus  in  place  of  (6,77)  s and  (6,75)  * we  have 
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* C,  Strachey  and  P„  J,  ’Wallis  - Hahn' 3 Functions  5^  (d)  and  IL,  (dO  - Phil,,  Wag. 

22>  87  (1946), 
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A 


obtained  respectively  from  (6.77)  and  (6.75)  by  using  the  definite  integrals 
f jo-  rf,J-C<£f2!!3!r  c*or>r/-f  ?(*&] 

?/?*>'*  <*  - f- ' -^[-i  ^ ' / m^7 


(6.111) 


[Note;  There  is  an  error  in  the  formula  of  the  reference  for  the  series 
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s cos<rp  <5 
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where  and  G(‘?-,^)  are  replaced  by  (6.11C)  above.  The  tern  in  G(,^,p) 

has  been  omitted. J • 


Putting  (6.110)  in  (6.102)  gives  the  general  new  form  for  SA  ('?)  , which  agrees  with 
Strachey  and  V.'allis, 
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(6.112) 


For  the  special  case  (6.103)  , (6.112)  becomes 
IT  d na 

Sp(a)z  -ajs/npj cotap dp  ■---  -J-Jf/-coS  —^arl  dp  j 


p = an  integer  (6,113) 


which  also  follows  directly  from  (6.104)  on  U35.ng  (6,110)  . 
Integration  of  (6.110)  as  in  (6.90)  and  (6.92)  gives 
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Putting  (6.114)  in  (6.3.07)  gives  a gensrai  formula  for  (</)  which  agrees  with  that 
of  Strachey  and  'Wallis,  namely 


Upto)  - 'dEB  sin  2 rrp[/s>  2 s*n7r a -rafted  ^7 1 Wfe/tCg 

_ds*ni5°  ./?  tfYgJ  /-  & - fp-^sinzp^-rr)  -7T S ir>2.  rrp] ccfa<t> dj 

4-s  L f-  >/n  pj  /0 


(6  11.5) 
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For  p integrals  but  P/o  non-integral-  (6*115)  yields  the  special  formula 


U/o)--S£  C rr-<p)s)n2pj  cotofi  c/^j 

*-  So 


p>  an  integer  (6„ll6) 


which  is  equivalent  to  (6„108)  * When  both  and  p are  integers,  finite  contribu- 

tions come  from  the  terms  in  (6.115)  containing  (s/n 27rp  co^~^)  and  ('■  '">  *77% *>■(”#,/)  , 
to  give 
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p and  Sy#  integral  (6»117) 


equivalent  to  (6*109)  * 
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THE  STATIC  PROBLEM  OF  THE  CHANGE  OF  CROSS  SECTION' 


(7a)  Formulation  and  Solution  of  the  Problem 

As  the  frequency  approaches  zero,  or  the  wave  length  gets  long,  the  solution 
to  the  temperature  wave  reflection  problem  at  the  junctj  m of  two  circular  rods  goes  over  into 
the  static  problem  of  steady  heat  flow  past  such  a junction  produced  by  a steady  temperature 
difference « In  fact,  the  latter  problem  can  bo  exactly  solved  by  use  of  the  equivalent 
circuit  for  the  junction,  although  the  picture  of  incident  and  reflected  waves  on  a 
transmission  line  has  broken  down. 


ff 


i 

- ! 


T,  \ rod  I 


rod  JT 


'T 


J-u 


Consider  the  static  problem  in  the 
figure,  where  the  cross  sections  at 
e«~L^,L2  are  sufficiently  far  from 
the  junction  plane  at  z-O,  that  the 


higher  static  modes , of  the  form 


Jj&r)e'rn,*/ 


will  be  completely  'lamped  out 


on  those  cross  sections.  The  temperature  at  ard  away  from  the  junction,  thus 

takes  the  form  T-'Az-B  in  I f'  -A’z-'B 
pal  static  mode,  and  aacisfies  Laplace'1 


takes  the  form  T-'A^  D in  1 zn*  A"  -A5  in  II,  which  it  the  form  of  the  first  or  princi- 


irten  we  must  nave 


■'  / ~ 
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while 


At  zo  , the  principal  static  mode  takes  the  form  TY°)z  ^ 

-4  JT'A  z a' 

^3-  ~ ^ holds  for  all  z including  z*»0.  Although  this  static 

mode  does  not  satisfy  transmission  line  equations,  we  may  consider  that  T and 
define  a voltage  and  a current  at  the  junction  plane  by  means  of  the  equations 
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•*n  using  = 'L/KnR‘-)t  ^ so  that  ^ *?<>  i independent  of  frequency. 

Similarly  Va  'T‘/oJy  Z'(o) -K */?'<(*' r/^)a  . How  the  four  quantities 
will  he  related  by  means  of  the  impedance  matrix  used  in  the  general  treatment  of  the 
junction,  and  in  fact  the  equivalent  circuit  may  be  employed,  but  with  the  static  forms 
of  the  circuit  elements,  i.  e.,  their  limits  as  oJ  o when  wa  shall  put  7 2s  , 


The  circuit  gives  -4  *■  Tc  since  it  consists  only  of  a series  element  or 
Krr#2J  * 

and  a.s  in  the  wave  problems 

or  ZO-  & ■'  Z-S 
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Thus  we  have  four  equations,  which  determine  A,  A/ & ' in  terms  of  the  bo  undary 

condition  and  the  Junction  impedance  2s  , hence  solve  the  static  problem. 
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The  temperature  varies  linearly  with  z from  T-.  at  z—Lj*  until  the  higher  inodes 
around  the  junction  have  significant  amplitudes*  the  principal  mode  gart  suffers  a 
discontinuity  at  the  junction,  then  varies  linearly  to  T?  at  z-Lg*  inc  hoat  flow  is 
given  by  -7?X  #*■/!  which  is  constant  over  the  entire  length  or  the  rod.  Since  this 
remains  the  same  even  near  the  junction,  ii  is  evident  that  the  higher  modes  make  no 
net  contribution  to  the  heat  flow.  This  result  follows  just  as  in  the  dynamic  problem 


immediately  from  the  relation 
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where  the  integral  is  proportional  to  the  heao  flow  through  a cross  section  due  to  the 
nth  higher  mode.  If  in  (7.5)  we  take  X - A*  / /?»/?'✓  2s  T° 

the  heat  flow  becomes  simply  -rr  X/ f*2  (Tx’r71)/A/^zzj  y the  result  for 

a uniform  rod. 

(7b)  Similarity  of  the  Static  and  Dynamic  Fields  in  the  Junction  Plane. 

The  basic  reason  why  it  has  been  possible  oo  analyze  the  static  problem  in  the 
same  terms  as  the  dynamic  problem,  is  that  the  situation  in  the  junction  plane  is  the 
same  for  both.  The  modes  of  both  problems  have  the  same  radial  structure,  and  for  a 
matching  problem  at  fixed  z (i \ junction  region  with  no  extension  in  the  z direction) 
it  dees  not  master  what  the  z dependence  is.  Thus  the  values  of  T and 
in  I and  II  arising  at  z«0  from  the  principal  static  mode  might-  be  thought  of  as  coming 
from  the  principal  modes  of  some  dynamic  problem  at  finite  frequency  ( which  have 
cc;*}  , dependence),  and  determine  the  same  higher  mode  structure  in  the 

junction  plane,  characterissd  by  2?j. 

One  could,  in  fact,  repeat  the  analysis  of  Sacticn  5 exactly  for  the  static  1 

problem,  using  only  di  'ferenfc  fortar-  for  the  principal,  mode  functions,  and  arrive  fit 
the  same  formulation.  This  is  namely  that  the  principal  mode  part  of 

is  continuous  across  the  junction,  and  could  identified  a3  the  current  ( it  is  the  total 
heat  flux),  and  that  the  principal  mode  part  of  T is  discontinuous  across  the  junction,  and 
on  each  side  could  ( bolder. vi mc-c:  as  tne  voltage.  'This  discontinuity  in  \A>  is 
proportional  to  the'„  -Asa at  ci  y*  because  of  the  linearity  of  the  problem,  and  the 

proportionality  constaav  repiosv.ts  „• ■.jsonci.ally  the  contributions  of  the  higher  modes  to 
the  temperature  field  at  tie  junction  x'or  a particular  value  of  the  heat  flow,  as  compared 
to  unit  heat  flow.  (The  total  value  of  T is,  of  course,  continuous  at  z»0) 


/ 
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Hi*  higher  mode  amplitudes  and  the  complete  temperature  field  of  the  static 
problem  could  be  found  by  solving  e linear  equations  for  the  An  by  Borne 
successive  approximation  method. 

(7c)  The  Static  Problem  with  a Free  End. 

Y/e  note  that  if  one  end  of  the  rod  should  be  free,  instead  of  being  at  a fixed 
temperature,  than  A' - * r/~ > * c.  Hence  by  (7*3)  A as  well, 

and  by  (7.4)  B*  . Thus  the  temperature  is  constant  throughout  the  rods, 
regardless  of  the  presence  of  the  junction,  and  this  problem  has  a trivial  solution. 
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3,  THE  RING  SOURCE  AND  ITS  CIRCUIT  RaPRSSENTA? ICN 
( 3a ) Description  of  the  Source  and  Green's  Theorem 

A more  easily  realized  source  than  the  uniform  discs  over  the  cross  section  used 
in  Section  3 is  a thin  ring  around  the  cylinder,  extending  a distance  £ along  the  cylinder, 
which  produce,*  heat  uniformly  over  its  area,.  Such  a source  is  closely  approximated  by  a 
thin  resistive  layer  spread  on  the  surface  of  the  cylinder,  through  which  electric  current 
passeso  The  current  ray  be  steady  or  have  any  frequency0  In  the  latter  case  the  heat 
produced  is  the  superposition  of  a steady  source  and  a sinusoidally  alternating  one  (going 
ailemauwiy  positive  and  negative)  as  shown  by  (8.21)  below. 

The  ring  source  will  generate  higher  modes  as  well  as  the  principal  mode  , but  these 
will  damp  out  in  a short  distance  leaving  only  the  principal  mode*  It  i»  of  interest  to 
calculate  the  magnitude  and  phase  of  this  principal  mode  part  for  a given  strength  ring 
source  and  compare  it  with  the  idealized  disc  type  source.  Also  the  range  and  strength 
of  the  next  higher  modes  are  useful  to  know  in  order  to  estimate  the  isolation  distance 
of  the  source,  beyond  which  it  may  be  treated  simply  as  a generator  of  principal  mode 
waves. 


The  temperature  field  produced  by  a 

steady  oscillating  ring  source  Q& 

(so  that  ( is  the  heat 

produced  per  unit  area  of  the  ring) is 

easily  found  from  the  appropriate  Green's  function  G(0  r'J  for  a heat  source  in  a 
circular  cylinder.  Applying  Creep's  theorem  to  the  functions  7~ff) and  over  the  entire 

volume  of  the  cylinder,  where  GCfyr'J  will  be  defined  explicitly  below,  gives 

J ([€(%? )(V%+*x)TfrJ  - T(f){r V <s ( r'falA  y^Jo'S  , (SJ) 


Here  n is  the  outward  normal,  V and  S the  volume  and  surface  area  of  the  cylinder,  V 1 
the  Lapl-aci&n  in  the  coordinates  r , and  the  infinite  extension  in  the  z direction 
causes  no  trouble  because  all  functions  are  exponentially  damped  away  from  the  source. 


(3b)  Determination  of  the  Green's  junction. 


Now  choose  Qj\tr*r~  J so  that  it  satisfies  the  wave  equation  in  f everywhere 
except  at  ? * r'  , where  it  satisfies  a source  equation.  ( The  primed  variables 
are  not  to  be  aenfuoed  with  their  previous  use  for  the  second  rod.  Only  one  rod  is  con- 
sidered hers.) 


Since  we  shall  be  interested  only  in  functions  independent  of  azimuth,  i.e., 
circularly  symmetric  about  the  z axis,  x*e  choose  the  source  to  be  a ring  source  concern 
trated  around  a ring  at  particular  values  of  r*  and  a',  of  such  strength  that  integrating 
it  over  the  entire  volume  gives  unity.  Thus,  for  the  source  at  3 let SfCJ-jr^g  ) 

be  the  value  at  r.z,  and 


(8.1) 
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Thus  S(rj^,r't  ')  is  a d function  which  may  bo  written 

Sdst.r-J)* 


where  is  the  ordinary  Dirac  delta  function. 


(S3) 


Now  choose  G(fj/'*)  to  satisfy 


i'wid  the  boundary  condition 
^ * 


-o 


Then,  noting  that  ( V 2 +&?*)  7~~  Q/  (8,1)  gives 

oa 
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Introduce  the  relation  between  the  source  ana  C^%r)  /?  namely 
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Thus  j_s  known  over  the  entire  surface,  which  is  the  reason  for  the  choice 

of  boundary  condition  (8. 5)  on  G,  and  (8.6)  gives  directly  a formula  for  the  ♦eraperature 
field  due  to  the  ring  source j 

* 


7 ^ / srn.is^rjdj.  . 
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Thus  the  problem  la  solved  when  the  Green* r>  "unction  is  known. 

New  <*{fs£y  ')  is  easily  found  as  & mods  8Xp4n^^-0ri|  if  v/q  note  thu u 

it  satisfies  the  ordinary  wave  equation  except  at  r»rf,  z*z* • Hence  wj  put 


Gfarij  9-ZJ,  (f,r)eiJC"l} '*  Un  <rS3 ') 
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where  the  z dependence  is  chosen  so  that  only  outgoing  waves  are  produced.  (8,10) 
satisfies  (8.4)  term  by  tern,  except  at  z*=z% . However,  if  (3*4)  is  integrated  over 
an  infinitesimal  range  of  z about  z',  we  obtain  an  expression  for  the  discontinuity 
of  (**/*})  at 3 "*■3- ' , cn  assumi ng  that  G and  it3  derivatives  with  respect  to  r az*e 
continuous  at  a*z'0  Thus 
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Putting  (8.10)  in  (3.11)  gives 
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Multiplying  by  rJo  (<//,  A'  )dr  and  integrating  from  0 to  R gives 

fn  (r‘*2-  'J  T-Jo  (*n  r'J//ZTryC& 
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Thus  we  see  that  -fn  ( J is  really  independent  of  a5.  Putting  it  in  (8.10)  gives 
the  desired  result  for  0,  namely 
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(Sc)  The  Explicit  Form  of  the  Tarrmeratnrv:i  Field  and  Ocmparison  with  Field  of  Disc  Source 

Now  that  (8.14)  gives  the  Greer's  function,  we  can  find  77//^  'J  by  putting  it 
in  (8.9)  to  obtain 
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Thus  the  magnitude  of  the  outgoing  lowest  mode  wave  (and  all  other  modes  } beyond  the 
source  region  is  explicitly  dotartr-inedo 
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Beyond  the  edge  of  the  source,  i„e,- , forl^V^-^  , the  principal  msxie  consists 
of*  an  outward  propegai-iog  «»vp  r>r  ruum-i t.dH« 
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and  the  corresponding  current  of  r.he  equivalent  transmission  line  is 
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Thus  w®  see  that  a fraction  •r/"l  *&)  of  the  heat  of  the  source 

goes  into  the  principal  mode  waves  propagating  in  the  two  directions  away  from  -he 
source,,  Ihis  fraction,  which  is  a typical  interference  factor,  is  very  near  unity 
when  )££<*•/  i.e.,  for  long  wavelengths  compared  to  the  extent  of  the  ring, 

as  will  be  true  for  most  ractical  cases. 


For  a simple  disc  source  of  c,he  kind  considered  in  Section  3,  generating  heat 
Q ergs  per  second  uniformly  over  the  cro^s  section,  only  principal  mode  waves  are  set 
up,  and  the  magnitude  is  given  by 
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This  differs  only  by  the  factor  sinf7£4*  )/'CjCJt/i)  in  phase  ant'  amplitude 

from  the  ring  source.  For  all  practical  purposes,  the  ring  sourco  is  thus  equivalent 
to  the  disc  source  for  production  of  the  temperature  wave.  However  there  are  also 
higher  mode  waves  generated  by  the  foricar,  and  it  is  desirable  to  get  beyond  the 
range  of  these c The  ratio  of  the  amplitude  of  the  first  higher  mode  wave  in  (3015) 
compared  to  the  principal  mode  is 

MM  JL  'x>*  s ¥ M - -7?  (8.  M 
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where  the  approximate  value  is  obtained  on  putting  r8®R,  X-  X,  and 

neglecting  &i  and  . Thus  , the  rat  io  is  ./z/?yj^2r  .0/  i~o  .00/ 


for  reasonable  values  of  R and  A 


distance  of  one  radius  from  the  source. 


and  thus  can  be  safely  neglected  at  a 


Finally  we  note  that  the  ratio  of  the  temps raiure  (principal  mode  part)  at 
position  s’  to  the  Instaneous  heat  flux  at  the  source  is 
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where  IX/  - f /&  , and  the  approximation  lXl£  l has 

been  used.  This  gives  the  relative  phase  of  the  temperature  at  a point  ?,*  compared 
with  the  phase  of  the  source  ( which  is  the  phase  of  the  second  harmonic  of  the 


electric  current  producing,  oho  heat  of  the  source  }„  The  pi  incipai  tern  close  to  the 


source  is  the  45°  lead  of  the  temperature  . At  £ ' - o 
impedance  of  the  infinite  rod  to  the  ring  source, 


it  gives  the  thermal  input 


"V 

> 


4* 


i w i 


4 

(tfd)  ieid  of  tbo  Si.at.ic  Ring  Scarce 


i'be  above  analysis  using  Lhe  Gr«.en*s  function  method  may  also  be  a;  plied 
solve  lhe  problem  of  uho  sialic  ring  source , ihis  is  of  interest  since  the  method 
of  electrical  heating  prouuces  such  a static  source  3uperposed  on  the  dynat.de  b? 
{*riodi>:  source^  Thus,  if  the  electric  current,  io  given  by 
,f<'edin r.  into  a resistance  /t!  . the  instantaneous  heat  produced  is 

= ^(,+  ces2u,.t)^  Q,  + a Q.X*'**  fSZO 


where  is  a steady  heat  source  and  Q_  the  amplitude  of  the  periodic  heat  source 
with  frequency  double  the  driving  electric  current.. 


fhe  static  Green’3  .function  will  be  defined  as  a solution  of  the 

inhomogeneous  Laplace  equation,  in  place  of  (8 .4), 


Fzc5 


j-(r-r') 

2.77  f 


so  V v 


(822) 


with  boundary  condition,  as  in  ( S .. 5 ) 
^6“ 


(8.23) 


Jr> 


- o 


Then,  just  as  in  < >’  ?o)  - (8.14),  but  using  static  mod»3  in  the  expansion 

Gy  ( r' 3 j O 3 ) ~ 8 / /e2  ~3  ''8(2',  3 ')  -+)>  1^4  (fr> •)&  ( /'j j ’) 

(8.24-) 

-OP  7. 


/ r 8 (Xnr.UJSn  r‘)  ^8,/ 3 3 fl 

‘ 7 *£  -£j7T^kJ—  e J 


f 

l 

L 


where  the  constant  C is  arbitrary.  Note  that  the  linear  term  in  % follows  from  the 
lowest  mode  term  of  (8d4)  in  the  1 bruit  X — * & 

how  let  the  lowest  mode  yurt  of  T(  X J-)  i or  /3 / , hence 

outside  the  source,  be  A a-/*Bv  as  in  section  7 and  consider  a length  2L  of  rod, 
the  ends  maintained  at  temperature  7,  3r.d  the  middle  having  the  steady  ring  source 
Ql  o Then  k and  S are  determined  by 


///.-/  3 = /, 
— = -t<rr/?\4 

P 


(3.25) 


x 

• ¥ 


(8  26) 


1 
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where  (3,2 b)  fo . lows  because  only  the  principal  mod©  carries  head  across  any  cross 
section. 


Now  apply  (301)  to  Gg  and  T within  the  rod  of  length  ?1.„  to  give 

rj 

r(r& j - f& (*.* 'irXgjl ”■*'* 


,/f 


i-fi  -JzC^fvrdr  -yf*  (- 3g)  2 rrrdr 


'C-L 


Putting 

(31 ) 9 Sj f0r  <a4  A ; ( ff) 

( Ztr/HK  z XL 

and  using  only  tho  lowest  mode  parts  of  G and  ( ~0/)  0?  ^ L gives 

rfr,„,  r%.  . -t-  r vsg'-A>;  * 

71 rnj,  /*.*)  e 'Mfie ’ 


~ a _f£Tj 

2 A'TT/?1  ‘ '/r  .f 


+jftc~  fzhr>-)Z7Trdr  +X*(C  ~ *££&)(’ ~i^c^)z'rrdr 


fa  ?#) 


-/  T,(-fFRi)  -f  77(-frfi)2’TI'dr  . j‘7-%  . 

Hence  the  solution  to  the  static  problem  is 

7W  ?>{  . r**>j 


Evidently  (Q„29)  has  the  correct  linear  behavior  in  the  lowest  mode,  and  agrees 
with  (3„15  a)  when  Kr>  trere  is  replaced  by  si  ft,  (Note  that  the 

arbitrary  constant  C has  cancelled  out  of  the  result,) 


V 


Appendix  I ~ LI3T  AND  DISFINITIOHS  OF  SYMBOLS  IN 
ORDER  OF  FIRST  APPEARANCE 


1)  Q 

2)  T(r,t) 

3)  K 

4)  f 

5)  Cg 

6)  Z? 

7)  cj 

8) 

9)  Tfo)(F) 

10)  X . k 

11)  %n 


13)  Jc(rnrJ, 
' U)  jSn 
15)  /? 

' 16)  77,  77 

17)  Qi 
is)  A* 

' 19)^^,  Zz 


20)  <3* 


m 


21)  7-7  7- 


(Cross  references  are  given  for  cases 
of  repeated  or  nearly  identical  symbols) 

vector  heat  flux  (jol"^c/tizs^c  ) (cf.  /J20, 4?  VO  ) 

temperature  field  as  a function  of  vector  position  and  time 
(degrees  centigra &s)(cf.  £ /&2£  46  ) 

thermal  conductivity  ( JOJ^s/cmsec  tfcj  ) 

density  ( 3 ) 

specific  heat  ( ^0<J  a1  eg  ) ( cf.  36,  J~g  s~s~) 

diffusivity  =7  ^//>  Cgj(C  ^/scc ) (cf 
angular  frequency  w 2 TTV 

cylindrical  radial  and  axial  coordinates  (c-f.  6oJ 


-X 

w . * (r)e 

(the  superscript  is  dropped  in  later  sett  Ions)  (cf.  2) 

principal  mode  propagation  constant-  =(pi)( ~(/i/.)k  where  A is  the 

, . real  watre  number  (cf/// 

nth  mode  (complex)  propagation  constant  in  factor 

<r  ^ **  y , jtJ  - r»  - * / *) 

a (real)  radial ”w»ve  number11  determined  by  the  condition  ^(t//,/r )*0 

fc  f.  sh  e/j 

Bessel  functions  of  order  0 and  1 of  radial  coordinate 

nth  root  of  u)  l?-)  * C>s  /?- <%  *,2j  ***  j hi  if  mA/i  (C  -f.  3 2,  8 fj 

radius  cf  circular  cylindi  luttl  i Od.  fc /.  2 S/J 

cteady  temperatures  used  as  boundary  conditions  (for  static  heat 
flux  problems) /c/.  2) 

magnitude  of  steady  heat  flux  uuo  to  a static  source  ( cf • ( 2 oj 
cross  section  area  of  general  cylindrical  rod  (c-f.  3 / ff) 
lengths  of  rod  in  heat,  f 3 -Cv/  pro  bjL^ii'is  i’c-f.  £9j 
amplitude  of  harmonically  varying  source  ^ 


static  and  dy-ayulo  tempera' tures  9 solutions  of  static  and 

dyo"j?=3  c heat  .f’c'j  oroblemr.  respect;.*  rl.yj  total  T* » ~~ea.^. 

(cf  2.9) 


-Ali/C 
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22) 

?3) 

1 

1 

24) 

R 

25) 

2> 

1 

26) 

S„  -xurt 

P** 

27) 

28) 

\r0e'lujk 

29) 

an)  — $ 

3D  A,  3 

I 

32) 

33)  Vfr),T(£) 


representation  of  rod  terminating  in  free  sndj  the, surfaces 
are  all  solid  lines 

representation  of  rod  with  constant  temperature  end  at  static 

value  72 

re  presentation  of  steady  heat  source  Q,  , disc  type 

representation  of  dynamic  boundary  temperature  zero  on  end  of  rod 
representation  of  dynamic  disc  type  source  Q0  & 
representation  ox'  combined  static  and  dynamic  source 


representation  of  oscillating  temperature  bcnindary 

condition 

representation  of  combine a steady  plus  oscillating  temperature 
l undary  condi  tion 

representation  of  serai- infinite  extension  of  rod  to  the  right 

amplitudes  of  travelling  principal  mode  waves  in  on**  ~J>. 
directions  respectively  ('cf.'d>¥-¥-yS~<fJ 

amplitudes  of  cosine  and  sine  standing  principal  mode 
temperature  waves  (cf.  /4y  9/) 


(complex)  voltage  and  current  cn  equivalent  transmission  line  (cf.34  39,  4oJ 


34)  7~o 

‘ 35)  2°,  K 


values  of  VfoJ  ai  d Z'oJ  t hence  usually  refer  to  reference  plane  (c-f.  33,  4o) 

characteristic  impedance  and  admittance  of  equivalent  transmission 
line  (ct.  3 J?  4/y  ¥-2y  *3 


36)  C constant  factor  introduced  in  defining  V($)  - C7~ / 

* later  taken  as  unity  fc-f,  S~y  S~oy  s~s) 

37)  complex  voltage  ( or  temperature)  reflection  coefficient 

at  a termination  fci.  /£  73) 

3&)Aj3j<*ij8i  Yfjlt,  quantities  referring  to  rod  IIS  as  opposed  to  rod  I,  in  genen-l 

the  prime  will  be  used  for  this  purpose  Ccf.  3fy  3Zy33,3¥,  33^36,  etc.) 


f- 

V 


39)Z„,22,2^,Zlt 
' W)  K ",  I." 

V 

f 


elements  of  impedance  matrix  relating  Z>y[£yl0'  (cf.  33y  4(42,  43) 

voltage  and  curi  vnt  defined  in  rod  II  with  opposite  sign  convention 

for  c-,rronta*  K "=  K> /„  "=  -/  •)  (cf.  3 ? 34) 


■I 


42) 


43) 


44)  Aft 

45)  £M 

46)  5&u*%i 

47)  <?c*7 


48)  -? 


*♦// 


50)  c£. 


m 


51) 

52)  =* 


53)  ^5), 

54)  W 
P&h  ¥"<}) 

55)  Ck 

56)  Ai 

57)  ^v? 

58)  4 JU 


59)  ^ £» 


60)  3y 
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.".pedaneo  raalrix  ej.enents  relating  K.  HZJL”  fc-f-  39J 


arbitrary  Impedance  terriinatlng  red  II  at  distance  X 
free  reference  p>  ax»  fcf-3%  *-3,  +3) 


impedance  of  rod  U seen  at  reference  place  = ( ’)  fc-f  3%  fz,  &3J 


amplitudes  of  higher  modes  in  rods  I and  13,  /?•  'Z,~ fcf.  /^3/J 
the  function  K(  'r/^^  ) in  the  Junction  plane  fs-  -a)  fe-f.  S’T’J 
principal  naode  parts  of  7”  and  (^^j)  fcf.Z) 


total  iconspisaihost-  flux  through  the  cross  section  of  the  rod 
at  position  fc-f./) 


series  ipodancc  rs presaati ng  the  effect-  of  the  chsngo  of 
cr???  section  on  transmission  lines  I and  II  fcf.  3^  4Z.+3J 


a BiLiErical  constant  depsnding  on  ( ^ >?'  is,  fn'  £?ut 
independent  of  frequency  ta.£) 


• nesaerlcsl  constant  demanding  on  sot;  4».i  -for**//*!  oot/onj^t 

(hence  on  frequency)  ( cf.  S~,  36,  S~S) 


the  ratio  /^*/  2/  /y  contse/tr/on.  fc~f.  &*) 

is  equal  to  } 2f  is  asymptotically  ib 


is  identically  equal  to  s ie  graatear  (less;  than  or  equal  to 

is  approfxiraateiy  equal  to  j C~)  is  fonaally  equal  to 

s*lsn,  (<**%** 

*93-9  ■** 


4UO  J>*3  function  or  the  logjiritieic  derivative  of  the  rectorial  ftoctica 
cter-*vathf«s  ;f  x7 ( 

Siler1  s constant  - - S~77 2 ff  ( c f.  S~,  36 J 


the  Bernoulli  smbers  fc~f-  3/) 
the  2uier  nabe rs  fcf  ^3" 

jl’&.X,  1.--  -tf/-  (‘S&'Y1  re*  mo 

s - Jr  Sn  S — X^JLT  ,>J  rs  c symjoXzf^fc  /«  Ss>  /or-  Xo/y*  /?  fcf-  gSj 


. I 


y=  y^-^=  (”%>)( i-r) 

not&tlc'T  rased  for  chwaIkk*  in  the  rsyn  fonr*aj*t"  cf  ssctlon  6 fcf  3) 


- 

v 

r% ... 


n 

1 


a *3 


,1 
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JZ""' 

61)  , *4. 

62)  £,0j£ 

63)  Sfa^l/^fa^) 

64)  O(j).  Sit}) 

65}  f**  0V  Jfi tw 

66)  xv 

67)  Gzfa}),  fz  la,}) 


gfJLvnl 

■fihus  *£•  ■=■ 

nsimericiu.  constants 


rTsP-,«l  JP  ^ jr  f&t~J 

J ^ 


-r“.  ^ 

frf  S3,SSJ 


- o- - ^ Jfe,  fzf  ^) 


functions  of  t»  variables  defined  by  snas  in  (6.65),  (6.66),  (6.69) 
and  required  in  asymptotic  form*  for  C*r  id  67) 

cosine  and  sine  integral  functions  ’ X ^ 

t"-**  ££*" - s„'e* 


asymptotic  values  of  //  a, 
£ « (<rf  62) 


jsr'jzy 

obtained  by  pu  ting  j9,  * '*■ 


functions  of  two  variables  defined  by  sums  in  6.89  ;«jjis8d 
in  asymptotic  fens  fear  Cy^  *•*) 


68)  S^foj , Cfafa)  tire  functions  defined  by  W,  Haim  as  courier  smo.givin  in  fc  f- S?) 


69)  i 

70)  <2 

?i)  Gfr  r') 


nftSCn^K  3 ■) 


73)  J(r-r-) 
1i^6£ff.P‘) 

75)  yf<5 

76)  £-n  (*) 
Tf  )£**>*  (*) 

78) 


length  of  ring  source 


uft 


strength  of  ring  scarce,  producing  heat  4 *7,20} 

Crccr.*  c function  ir.  infinite  circular  cylinder  with  free  wall  boundary 
conditions  (^/?st}=o)  fc{-  ~7+) 

a ring  source  func  ten  giving  strength  at  ft j, l,  cf  a cf&/ /u 
function  ring  at  Kr } . Ccf.  33  *7e-J 

the  /?iroc  delta  function  fet  Btj 

the  static  Green's  function  corresponding  to  6/f//'y  fc-f-  7fJ 

real,  imaginary  rerts  of  , ^ 

tin  innsml  C/nOS'e***** 

the  incomplete  eaqpooantial  of  ^ . » consisting  of  the  first  /7 

tens  of  mm3  power  series  ii}/ 

/“Ti  ^ 


the  eoffipJeaent  cT  the  inemrstte  eacp^nential-  consisting  of  the  tail 

of  the  power  series  fnm  <m 


Bh 


i 


79)  [ fJ 

SO)  S*oh/'& 

CoS,  (TC&  Ssn*  /*} 


tfca  largest  integer  is  s7z 

iacanpiete  cosine  and  si^c  fuuiciioiH  and  ".h&JLr  c crap leases ts,  defined 
in  (II  .8), 


81)  T/^?/,  a rf  ISj  A 4 de  pj  ^ °0^  «CT  ^ -G 

are  all  sryrabols  occurring  with  special  meanings  In  Appendix  Xlb^ 
dif lerent  froa  Use  uses  in  tltei  i!iin  text  and  defined  above. 

For  these  special  definitions  refer  to  lib. 


Appendix  II  — EVALUATION  OF  INTEGRALS  IN  SUM  FQAKULAS  f 0-1  THE  Cj 

22. 


(II&)  Tcmbt  Series  Expansion  Method 

The  integrals  in  (6.69)  may  be  evaluated  by  expanding  cot 
in  powers  of  s and  integratiiig  term  by  term.  ihe  resulting 
expansion  is  convergent  cut  to  H'?71'  which  is,  in  fact, 
the  largest  value  of  £ - 2 7T//-  that  occurs  in  the  change  of 
cross  section  problem.  However,  the  expansion  need  never  be 
used  for  values  of  7,  greater  than  zr  f where  it  is  strongly 
convergent,  since  we  have  the  identity 

O,  . J.7T-J. 


Je~3^catd^  a 3-  - &AQScaf&.  <*3. 


(2Z./J 


Her  values  of  greater  than  . which  may  occur  in  other 
applications  (see  the  section  on  the  Hahn  functions)  the  periodicity 
of  cor3-/z  (period  Zrr  ) permits  reduction  of 

(z  /nHfrr  77*  (2/n-rZ)rr  Zrr 

J to  J and  <rf  J to  J 

ZmTT  frmujtr  7T  ^ 

which  by  (11,1)  may  also  be  reduced  to  Jo  . 

In  fact  3-  may  be  reduced  in  the  original  arguments  of  Gfcte) 

by  the  easily  derived  formulas,  where  . is  en  integer, 

/rfOjx)-  cos  2 rrop/^fa £ -Zrrp)-  s/r>2>< ..  ^GCa^-Z/rp) 

(2L.  /*) 

£(a>/ p*  cos  l Trap  -2 np)  P sin 2.  rrctp F(oy£  -Zrrp) 

It  is  also  noteworthy  that  va  ues  of  a greater  than  / need 
not  be  considered,  since  the  original  arguments  in 
may  be  reduced  by  means  of 

cos(cr-<jf)3~ 

P &,?)--  ?(«-*> 3 1 - Z.  ~ 

(JZ.  fc) 

6U})  --  <?(*-/»})  - jf 


l 


It  is  useful  to  reduce  3-  and  <7  in  this  way  since  the  con- 
vergence of  th«  series  (II o 10)  and  (11,12)  is  then  improved  by 
making  the  quantities  S/n^fO p and  Ccarfop.)  decrease  more  rapidly 
with  D o Also  the  integration  range  in  the  numerical  integration 
method  of  Aooendix  Tib  is  reduced. 


Thus  we  have  , 

jfeia*cet%  --  f £,  "(i/ft  (t)  )<Tp 

^ t 

: [Cl  (op  ~Q(c;t)Jp  J,(Sj(ap - S/,'arr;)-2 Je 


(JZ.  2) 


I 


9 


I 


The  ints^rais  in  (II.2  ) nay  be  explicitly  evaluated  In 
simple  series  fora  as  Xoiluwa 

— . ; * — fl,  LfcmH 

(2/)-# Je  S aznL  ^ J 


atmre 

4,'V*  A? {****”<*  = '*e * % 
Ijti*  -J.ee**-*) 


(zz?) 


Repeated  applies ties  of  (11.4)  yields 

. /J  ■ $ _ 

mi/-  ii  r-  ~y-  nr.  i 


/.  ah  * (/-&  4L  * — / - -t  // 

^/l  1 //  J 1 


fz-e1**  Jr* 


c-*)] 


(xr.s-j 


«h«r«  we  define  the  incomplete  sapeasst  fs"^  and  its 

coapienaent  <*-  Hy 

^ "T7  ^ ^X~ 

a*  / 


4zr.  <*7 


rt 


cw: 

-/E 


^ 2L4 

x/7«  j^/ 


For  the  imaginary  argument  in  (JI-5)  it  is  convenient  to 
separate  ^v-/v'  into  Ute  lacuaplet-?  cosine  and  sin©  functions^ 
defined  analogously  by 


JE*rP  ' 

srrt 


* atof,-*atvjw  t** 7) 


where  the  notation  if/  ssuss  the  largest  integer  in 
Thus  we  have  ^ 

G*  '*/  F J /v'A,  = — S,-/^2".  - 

**•  s*o  fl”  4,*n+i  fW-‘ 


/s*n*i  f ^/V-' 


;>  <*■  _ ~r  , , 

r.  fx)  - y f~)P — *sinX-—S  f~)P -S-mjt-Snrt 

jm~  -£  J te+os  fr„,.  <V.hV  - »« 


6z^ 


Putting  (II„5)  into  (II.3)  gives  for  the  Integrals  i Xk'ZLZ) 


e 
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X 


Finally  we  can  write  ihe  real  and  imaginary  parts  of  (IT  2^ 
which  occur  respectively  in  ST5C»  ^4/^  ? in  the  forms 

(II  9)  and  (11,11)  , 


(zz.fj 


fc~saJ'  Cttj/aj^tesjtS^tojj-crMrC^fi^-S'rforr  JT^, 
n-i  0**L 


Since  the  functions  Cst>  and  S*n^  do  not  become  1 
as  JO -*■<£>  but  approach  cos  and  tin  , respectively,  it 

is  best  to  rewrite  (II  9)  in  terms  tif  the  complementary  functions, 
which  go  rapidly  to  aero 


- C/faj)  —Csfan) 


( ZZ.  /&J 


(lib) 


Similarly  ^ /3u,f-)JW 

* Sj.'Oj )- &(<>»)  ~ sL,  &*/&/(&  4 

**  -s  ,nc>rCc^:'rrrJ+<rcx^J^**J/ 

t fjfiyz ) *’ ria?'  f**^*})*  Ccsm^^^fJ  -f-SifiorrChi,  fa^-d3ar£nrSjfftJm»rj 

x-  (JJL/2) 

Itefrical  Integration  Kethod  and  Fllon*  s Formula 


The  evaluation  of  the  integrals  in  (6  69)  or*  better,  of  the 
integral  in  (6,71) 


or 


= / e***  ff})*}- 


(ZT-'3J 


may  be  simply  and  accurately  dene  by  direct  nues-rical  integration, 
even  for  large  values  of  a 3 provided  we*  nake  use  of 

Filon*  s method  -^'whicn  treats  t-ie  trigonaiietric  factors  analytical !y , 
The  method  Is  based  on  Simpson's  rule,  which  requires  a tabulation 
of  the  integrand  fine  enough  so  that  the  integrand  nay  be  represented 
to  the  desired  accuracy  by  a parabola  through  each  consecutive 
triplet  of  points 

TFT 

' 'L  » G -lion.  Pro-  Pey  -e?  Sdin  hg  38  48  (I92S-9) 

See  also  C,  J Tranter  "integral  Transforms  in  Ma theoatlcqi  Physics?" 
pethsen  1951*  p 67- 
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Accordingly,  if  this  tabulation  is  done  for  ffj-)  , the  integral 

(IT-13)  may  be  evaluat'd  as  accurately  as  Simpson ‘ s rule  applied 
to  f(2-)  alone,  using  the  same  set  of  points,  by  doing 

analytically  tb  t integral  of  the  product  of  and  the  approxi- 

mating parabola-  The  analysis  is  repeated  here  briefly  for  completeness 
an  1 ease  c'-  re'erenc®,  and  slightly  extended 


Consider 

I = £ "e  * > *3  5 tk,2***' 


(J2.  ■'*) 


where  the  integration  range  £ to  d has  boon  divided  into  2/7 
each  of  width  /)  and  the  notation 

£-s  - c-f-sA  , r *Oj Z/> 


parts 
UZ.  Arj 


has  been  introduced  for  the  S**  abscissa.  The  integral  is  split  into 
H sub -integrals,  and  the  sub-integral  from  aue  to  2-nx-z  is  called  ■ 

Let  the  subdivision  of  (o^cj  be  fine  enough  so  tf.ai  in  each  of  the 
subintegrals  -f(j-)  is  represented  to  the  desired  accuracy  by  a parabola, 
hence  in  Xu*-/  we  may  use  the  three  point  lagianglan  formula  for  f(j)j 

•f( c y-  (stfi)A)  « rlU'lu  (1L./U 

*-  2. 


s-zi-h!  ; ~ / 4 r>  < / 


In  (11,16)  the  notation 
■fs  = f/>r>  f(c+S*) 


(2Z./7) 


has  been  used,  and  P is  the  usual  Lagrangian  coordinate  denoting  a 
running  fraction  of  the  basic  tabular  interval  /), 


Then 


T 

•*-s 


•t; 


(JZ.  /8J 


, e**tc+*Hh ) - (£i - ^ * (£s±*l  fst/ ? 


J 


2 


(z?:.uj 
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where* 


L e 


OP^/  = 2J£2± 


&•  #/) 


£,  * fe**/",,#  , 


F , f'e^e^e-^  gx/,g^  t£2t?-  ±i»2i  I 

2 “/  <5>  ^2 


~2a  ^->' <£*  , 2d  i/>7<9 


<5> 


<!P* 


(22.  Zo) 


Putting  (I~-PO)  into  (II. 19)  and  using  the  notation 

/J*  Mfr)  --  e'°*  ffr.)  a e^^pcfif,) 

for  the  complete  integrand  in(II  13)-  we  have 

^ 7 

-jr-+*jp*rj 


4*  e ^5  ■ ’ ~ 7 

4.  ** 


where 


oc  = 


_ Jt.->  <g/F4  • cos&  P,  •_  _4  s/n&cos&  _ 2 


P 


y-^ 

2 2 


«92. 


^ = ras£>4  Shi6/rt  * (/+CQS*&)  ~ 4-J1LL0SM& 

<02.  ^ ^ ?/ 


£7- 


<9  3 


.v--j  4-s'ne  _ jl  C&S& 

r-  « g2  -?r 


i 

y 


2 0 


(JT.22J 


(ZZ.  23) 


andety^  ^ are  the  3arae  quantities  used  by  Filers  or  Tranter 

Taking  the  sum  of  J~s  * . over  all  the  aubintegrals,  we 

have  for  the  entire  integral 

J,  %%„ -hfepfii  JfagKfcffiiw  *■** (a.  2 +) 


£=<. 


L . 


(lilt)  is  the  desired  integration  formula,  where  the  first  sum 
is  over  the  integrands  at  all  even  abscissas  minus  one  half  the 
first  and  last  i.  a * v;  (o^-  0 fc.  - 1 f 0 > * 1 for  &*0,  is  the 
Krone*  ek.  cr  delta  .symbol) , and  the  aecond  sum  is  over  the  inte- 
grands at  ail  odd  abscissas*, 

In  taking  the  real  and  imaginary  parts  of  (11-24)  to  evaluate 
the  integral  with  cosine  and  sine  factors,  the  last  term  with  the 
coefficient  jla  mixes  the  two  types  of  integrands  „ Thus  putting 


r.-  _ 


i-CV< 


r t ■, 

T ifc) 


J.j 


o*> 


C oltl  l,\ 


cfS 


* C, 


-V  i. 


Ur- 25) 


we  have 
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/7-y 

■ Afei' (Clt  - -«■(*  *-W ] Ctr.  2fJ 


pin ai  ffyjj).  - hJ}9jL  (szi~ 
C £ SO 


efe^cCrr  ^ r +r*  irL- 1“-  *>  fH,  Z7) 

— J£0  ~V*W  ' ~ • — ' 

2.  * £no 


Finally , we  note  that  if  & is  small,  (11,25)  is  inconvenient 
for  calculation  of  <f~  because  the  difference  of  large  nearly 
equal  terms  occurs  As  Tranter  indicates,  the  power  series  ex- 

pansions are  convenient  in  this  case,  and  we  may  write  these  in  the 
general  forms 


or  - 5 ' (~)n  

%lD  (lni-6)! 


(Z/l+2)  z.n+j 


_ ^ ....r,  (Zn-:)2z»^  &2„ 


/i-O 


(Z/>+3)t 


( 22. 2SJ 


f-  2 (-)n  ±l2p±iL  fit* 

(2n+3^{ 


/?  :o 


The  numerical  values  of  the  coefficients  are  tabulated  below  and  suffice 
in,  at  most,  ten  terms  to  calculate  to  6'  decimals 

At  most  six  terms  are  needed  for  For  6"?/  no  numerical  diffi- 

culties arise  in  using  (II,23)» 

In  applying  (11^24)  or  (11,26)  and  (II  27)  to  the  integral  (11-13^ 
we  have  merely  to  tabulate  ffy)*  co£3/z  at  small  enough 

intervals  over  the  range  of  interest  ~(07K*>  TTJ  , that  can  be  well 

represented  by  a parabola  over  these  intervals  Otherwise  stated, 
only  the  second  differences  of  the  tabulated  values  of  Pf*)  are  then 
needed  in  interpolating  Pf1})  within  the  range  As  shown  by  the  values 
helow,  if  we  choose  the  points  so  that  P *»  TT/Zo*  157030, 

then  the  third  and  fourth  differences  contribute  at  most  2 in  the  6th 
decimal,  and  that  only  near  gs-rr, 

rrr  r71  ^ 

Thus  the  procedure  is  this,,  Put,£  + Jj,  where 

is  the  hear- st.  even  multiple  of  tyf,  Then  if,  for  example, 

the  integral  contains  the  cosine,  as  in  (II  26),  tabulate  caja j*  Pf?s) 
betweer  and  7T  „ and  also  4 = S/n  f(3»)  a 

The  two  sum^  in  (II. 26)  containing,  at  most,  some  ten  terms  each, can  nuw 
be  added  up,  Finally  find  &-(a7r/%v)  and  c*,#,/'  and  substitute 
in  (It  26) o This  gives 
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y \ 

The  end  correction  JsoC<:tsa3'  t'st’t'-'Jb  requires  a second 

application  of  (11,26) » or  better  of  (l"t22)  since  only  one 
parabolic  5eci,iuu  Is  iWcJru  (>■  s 'v'io  so  that  is 

parabolic  over-  the  whole  range  -hog-.)  For  this  case  we 
hrv*  , e » a(t both  dif  ferent  from  the  main  integral,  and 

--  /r;_,  - //-;  9/t*  -F(*j&-) 

are  found  by  interpolation  from  the  table  of  , using 

second  differences. 


The  wrk  involved  in  tabulating  (II.  13)  by  direct  integration 
is  less,  in  generals  than  in  using  formulas  (II, 10)  and  (11*12) 
because  of  the  labor  in evaluating  the  sequences  of  the  four 
coefficients  / a , CosnfonJj  S/r>„/c/TrJ  „ 

Each  of  these  must  be  found  from  a separate  series*  whereas  the 
terms  in  (11,26)  require  only  tabulation  of  casof,  ( or sm 
cndntultiplication  by  the  known  values  of  ■f^)  „ ( Of  course 

^8,  must  also  be  evaluated  s and  the  end  correction  found) 
Ho£»ever?  if  more  decimal  d1*c«3  were  required  in  the  result* 
direct  integration  would  require  more  tabulation  of  -f fa)  and 

more  terms  in  the  sums,  so  that  eventually  the  series  (11,10)  and 
(11,12)  would  be  preferable  , Six  decimal,  places  should  be 
ample  for  any  practical  needs*  however. 


TABLE  OF  COEFFICIENTS  IN  EXPANSION  OF  <*,£ 


{a  } 


(8) 


UT) 


n 


(Zn+1)?. 


. ,A+,  li 

L-J  fzr:+3)! 


(■  ~ / (2  m-3)J 


0 

,044  144 

4 

„66b  666 

7 

1,333  3333 

1 

—.,006  349 

21 

1QO  Oil 

-/ 

-133  3333 

2 

,000  423 

280 

-c038  095 

24 

,004  76190 

3 

-,000  017 

102 

003  527 

337 

-oOCC  083183 

4 

.000  000 

/Aoo 

H'*  / w 

— . 000  179 

573 

001002 

5 

-,000  000 

00940 

,000  005 

9200 

—,000  000004 7 

6 

,000  000 

00014 

-,.000  000 

1378 

000  00000004 

7 

- ,000  QOO 

000  002 

000  000 

002  40 

8 

- 000  000 

000  0323 

9 

,000  000 

000  000  35 

FOR  NUMERICAL  INTEGRATION 


TART.F.  OF  VALUES  OF  £e^\ 

0F  f™*4,0*  2-& rr 


It 

S 

A?"  ~ 

n 

_ ATT 

3*  2p 

i~i  co+%. 

a'-'-./g+a" 

0 

0 

0 

0 

,157  080 

,013  095 

33 

2 

314  159 

,026  223 

65 

3 

c471  239 

,.039  U6 

97 

4 

628  319 

,052  707 

135 

5 

785  398 

,066  133 

168 

6 

.942  478 

079  728 

203 

7 

1 :099  558 

.093  532 

247 

s 

1 256  637 

107  584 

293 

A 

7 

1 .413  717 

> *S«» 

t 

343 

10 

1 570  796 

136  620 

394 

1J 

1 .727  876 

151  705 

454 

12 

1 ,884  956 

,167  245 

521 

13 

2 ,042  035 

183  308 

592 

1 i 

2 ,199  115 

199  965 

681 

15 

2 356  195 

217  306 

778 

16 

2 ,513  274 

,235  423 

88  9 

17 

2 670  354 

,254  442 

1C26 

18 

2 827  434 

,274  486 

1176 

19 

2 o9S4  513 

c295  712 

1365 

20 

3 *141  593 

c318  310 

1586 

Interpolation  formula  (Everett) 

g.fl*  +pA)  - <7  +£•'  +JF,  *vt£  " 

A"*-i  qfix-i)  - g.o*§.C|o) 

£<"  ~~T  « §-'*3-(^«+h) 

Use  of  this  formula  ^<ith  ths  tabulated  modified  second  differences  /V  gives 
six  decimals  in  over  entire  range  of  o£  5L 4rr  . 

Neglect  of  third  differences  makes  a maximum  error  of  2 in  ths  sixth  place; 
neglect  of  fourth  differences  a maximum  error  of  1 in  the  sixth  place.. 
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Appendix  HI 


Numerical  example  of  summation  of  asymptotic  form  of  the  higher 
order  ter® 


F, rom  (6-6o)  , x 

SOr)//>+Mj-J 

%>V'0-J£T\ 


f / “* — 1 

I (nt‘A)z} 


w/prr#  /?  - ~~ j **  ' 


Hence  by  (S.83) 


5- w ~r~\ \[~ ~ ¥(£)+-  &fi) + £ 4>  - £J-S  ^ 2 <^4  4/ 


where  b j (6.  6 S') 

g(u)- 


c>s£rr[-)-  #(*£)  + £ &(%r)]-£2[cosZ3.  rco.brrj- jJcasZj  corl 


and  by  (Z2.  /a) 


1 -V  „ „ , -V.  . °?  /'_)/V?_  f y s~  . _ . / (7  /r-0  > 

££rcoS^3'co^~-  ^3-*  Of 1 "C>i  ^ /-fc-w ^ caJ/i ' J 

PC  OS  kTrCasjS77)  *-s/r>&77  J/>^  f injJ 

where  by  (ZL  ^ 

/?*'/  “ 

Sr/^, 


> 


___  c"  -V  £/• 

=•  ^ r-^  — - = <•«•* 

/>:n 


Ml 


f *///! 


(%)  =ZO)F; 


<*  p 


/?■/ 


/?=/?  (lprt)( 


- S / s>  X-  — C—J 


p 


P'0 


(2/?p;7f 


Taking  r-J,  /?•=  ^ » = *>  ^ ^ , \KjZl  1 ifgfjfjgD 


ft"  •"««]  V 


„ 2.  TT 


sum  formula  -for  Jrit>'^  are  boxed  in) 

*}ss-*i>fr*2-*f  -l  &»(¥/* 

(j(+,S'  Z i r 3 /»=/  2.,  ten)!  *1  3/ 

~ /.  T7079S  -4-*.rooooo  + .693/4-7 -‘V+6'Z/-.  C5772Z 

— . JZc/oo 


(teiinc  that  enter  directly  into  the 


2/? 


kr-z)B„ 

n Z»Un)i  \ 

"l  .033  333  1 096  622 


&(9hPt*)  W7hv>7g?  S&'-PO-vAiSSS. ^&$*EM3L 


j*& sIzSI 


y cos  c- 


cs  T-f  i 


fraOisrj ; 


,004  861  1 202  580 

000  342  1,318  7« 

000  026  1 446  23 

,000002  1,58598 

.ooc  0Q02  1.735  23 


f .1 
1 


1 4* 


fT?! 


V4 


/rfaij)  - /T^  2-/J  = costal- j rW+i?'fi)j-zL^  ecs^c-^dj. 


•vf 

-~  f cos  ri  cof£  d\  ~ C/  ( %r  / ~r  '(* 

l/rr  " Z ' 


,-}  . *“  C~)r'cJ''  sJJr~  ,'8JI)  - £.T  ~ /£,,  ! 

/r  S. 7r-  / ^*3  ^-V  z /-■>"’  J J'^1  ~y  / 

ff-i  2.  /I  -r- 

y.  L-^.  y,r  Cos#  '*  rrj  y jv/i  4-7T  ( 4-f/J  } 


{-fS/n„  ( y^)  (~)  n CoSrt  (4'Tr) 


1 ,52083  ( -2) 

2 ,32552  (-4) 

3 96875  (-6) 

4 63583  (-7) 

5 72250  (-8) 


+1.500  000 
3-359  194  (1) 
1716  490  (2) 
3,08504  (2) 

2 93264  (2) 

1 76007  (2) 


2JT 


,12572  (-3) 

‘t 

tpowar  o-f  ten ) 


/ / 1 . , , , 

' Cos  crj.  = . /oy+33 -(-• 


cos  ®JT  ~ ~.  S'°°e,°c‘ 


7,511557 

.,000  000 

9 048334 

0) 

7,895  635 

(2) 

Sin  QIZ  ■=  . 

2 53^00 

(/' 

9-60074 

J 

3.21247 

(?) 

6 42932 

(3) 

Cos4-fT  *■  / 

2 38904 

(7) 

8,99333 

(3) 

s/n4-n  -0 

3 30672 

(4) 

00 HZ./ ) <~t>oooo(.oo90  i>ry  ' - 3 *6o2 s(.  C4-2 3 3 S~J /~.  OOi/9c>  +■  O 


- , O&S/Z  Q 


F<'*/'})*  S 000c o (, 4Z Z 7i34-)-i-  Co 0000 C ~?03/.r7)-;2ScoO{-  S~oooooFj)-  o&3 /2 & ■/-.oaa‘^2/' 


Si'  * ’ *rly 

f(1}>  • Ffb  ’/>  "sfl-iPGHWffl-yl™  yi+ *•*?]- 

j cei  J/wrrfi  rf?  - a(Y)-S>( ¥)  + K T^C7j^tC‘,‘  ?<r°-”^'77-  »«  Y^YY) 

y ccj  in  Cos,. [*£)+.  sin  ?zr Si n„  ( ^ )J 

s . / of  $2o  —(—  C&srtf-J  - JsDocoo  f.  oojroS/4-J  36602.5-  C 0+7/T3 7)  fo  — - OS6  3/Z 


r . OS- SYS' 7 


fi 

- — /T.T1 

c-rccu„irJ 

c-fjis>j¥/ 

r rr  < 

?n/3JT)Xn 

Si/i » ( 

* * 5‘oooco 

1 

„ (>80272  ( -2) 

+,500  000 

^6" 4S436 

+1,199558 

u) 

2 

0 55532 

(-4) 

2 636  726(1) 

5 .913485  U) 

2 09 5708 

(2) 

r/^72T  _ .3  66&2.S 

•*  r t / — ^ — 

3 

,21586 

(-5) 

9 39410  (1) 

1,17196  (2) 

1 129326 

(3) 

4 

18505 

( -6) 

1 21549  (2) 

1 03465  (?) 

2 .725819 

(35 

8 

6 

27464 

,62419 

(•  7) 
(-3; 

8 5224  (1) 

3 3229  (1) 

f 7943  (? ) 

• 743534 

C 3 

' ^ 

J/7»  r -/ 

4> 

Hence-  

)-  0-./+2 esrf- ^ ^°ooa^oj - . of-s'ij' 7 - / 


Thus  finally 

.0*33  773  A,  03SS-SS-L~  '•  ^ 53307/ *-Z(~.2Z7*&J  +-.2SS~7t4 

— ,OOS’*-+2  * * !2CC>3<  —,5Oc>0t>°(.2QS~7/f-}Th2(,/2.o/c>ci)J 


* ~*°+&37?3X  ./0//7/  « - .00+B9+5 

Direct  summation  gives  for  9 terns  .&• 
r?  s^°  S„(o>0  "*f 


1 

-*,001  548 

00366 

2 

",004  778 

- .00541 

e, 

5 

002  127 

* .00184 

4 

,001  423 

,00160 

5 

-001  366 

00133 

6 

0O0  155 

,00013 

7 

,000  087 

00009 

8 

,000  217 

COO  21 

9 

,000  037 

,00004 

£j'  ~ on f 7 ■&  2.  This  is  smaller  than  the 

exact  sum*  as  expected*  since  the 
later  terms  are  all  positive  and 
still  changing  in  the  fourth  place. 
The  result  also  illustrates  the 
difficulty  with  direct  summation, 
since  there  is  still  an  error  of  -4-fZ 
after  using  9 terms,  and  the  con- 
v*rgence  is  getting  slower, 
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APRP\rD IX  I«  .ViwricaJ  Tables  Useful  in  the  Calculation  of  ^ 


Table  I,  (a)  Bessel  Function  Hoots  and  Useful  Functions  of  the  Roots 


/-‘r> 


? rr^fnt^)* 


0 

0 

** 

3-83171 

068111 

109580 

.0082564 

2 

7 01559 

0203 18 

-037153 

0014157 

3 

10  17347 

009662 

015231 

,000469? 

4 

13  32369 

005633 

003867 

0002101 

5 

16-47063 

003686 

■005798 

. 0001114 

6 

19  61586 

302599 

. 004086 

0000661 

7 

22,76008 

001930 

.003034 

0:100473 

a 

25  9036? 

001490 

.002343 

-0000287 

9 

29.04683 

001185 

0018 63 

0000204 

10 

32  .18963 

000965 

001516 

0000150 

11 

35,33231 

000801. 

00:. 2 59 

00001.13 

I? 

38  .47477 

000676 

Of  *061 

0000083 

13 

41  61709 

,000577 

Gv, 09072 

00000*69 

14 

44  7S932 

000499 

0007842 

OOOOO56 

15 

47  90147 

0006847 

16 

51,04354 

.0006030 

17 

54  ' 8555 

C005351 

18 

57-32752 

000 ',780 

19 

60  46946 

, 0004296 

20 

63  61.136 

.0003882 

21 

66-75323 

0003525 

22 

69-89507 

0003216 

23 

73,03690 

.0002945 

24 

76,17870 

« 0002707 

25 

7932049 

. 0002497 
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•4  1- 
< 

\ 


n 

0 

1 

2 

3 

4 

5 

6 
7 


APPENDIX  IV, 

Table  I*  (b)  CoeTficients  in  Series  Expansions  if 


(lzn-z)  So 
2.f>  (Z/ 77-2)/ 

(2nt3)/ 

, 166  666  667 

>006  944  444 

»008  333  333 

-000  162  C37 

oOOO  992  064 

„000  006  102 

,000  168  100 

.000  000  292 

0G0  034  697 

.000  000  016 

,000  008  113 

..000  000  001 

000  002  067 

.000  000  561 

* 


its  Asymptotic 


Table  II  (a) 


Values  of  Basic  Sura 
Series  $lo‘rt(IT)  for  Various 


and 

r 


f 

X 

0 

0 

lo2 

o 0123  5 5 

.011947 

1,4 

02463 

02288 

1,6 

,03269 

,02884 

2 ob 

,04018 

03155 

2,5 

,04193 

02809 

3 0 

,040622 

023184 

AO 

03594 

01566 

5 0 

,03146 

c 0123 8 

6-0 

, 02769 

01244 

8,0 

02222 

01904 

10.0 

,01845 

.03075 

Zfao)  ,r  (f-  s/n(n+'A)  Z/£-) 
^ 'P/73  ' (S-‘t  Z+J3 


DO 


. fa  o) 


/~?c 


s'nr)  » ^ 4T  '(r) 

/!*/ 


S'^W' 


,n-t 


; 


-<r> 
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Tabic  II  (b)  Terra  by  Term  values  of  Basic  Sums  <c {aoi(f)  5nd  S(c‘,>(j') 
with  partial  and  total  sums 


n 

J„  U 2) 

1 JR  /4'W 

T 1 

r'w 

1 

,007611 

.007313 

j 

.02027 

01877 

.02979 

02617 

: 

2 

.002972 

.002900 

AAO  i C* 

'j  **0 

.00322 

.00126 

o 00101 

3 

,001131 

001108 

00008 

0000? 

.00057 

.00061 

1 

000326 

,000317 

. 00022 

00023 

! 

i .00063 

! 

.00062 

1 

1 

5 

.000011 

.,000039 

,0003' 

. 00031 

.00001 

.00000 

6 

o 000002 

,000003 

,.00007 

00007 

.00016 

.00016 

n 

i 

.,000037 

.000038 

i 

0000.1 

00001 

.00008 

.00008 

8 

000059 

.000059 

,00006 

i 

00007 

.00001 

.00001 

i_ 

-000018 

000018 

00005 

00005 

.00006 

.00006 

J"  = oizns 

S.O/  fa  3 3 

4s 

J>. 01072 

i 

C - .eits&r 

V/  " 

S*  PUZt-Y 

Js.oziaa 

i 

.02304- 

if-: 


i : 


•'V  ;Vl 


Table  lUb)  (ccmt.) 


n 

j£  /Z3  J*  / *■) 

srhzs)  jZ'~hsi 

S/^f3) 

.y/r  ( 

1 

s<c*(+) 

1 

0 

i 

03695  02819 

-,03465 

0 02064 

,029284 

,,013385 

o01939 

0 00251 

0 01387  ,00000 

2 

o 00058  00083 

00534 

,00562 

.009131 

„ 003494 

,01083 

.00783 

,00946  * 00489 

3 

,00172  00160 

00013 

,oocc6 

,,000510 

,000705 

,,00360 

.00362 

,00502  ,.00425 

4 

,00010  00012 

00106 

,00102 

,000409 

,000315 

0 00038 

,,00052 

.00181  -00190 

5 

.00040  00038 

0 00010 

,00011 

. 000678 

000668 

00007 

0 00003 

00029  0 00038 

6 

.00003  00004 

0001 3 

/v^^i  n 

yuux  f 

,000099 

r 00041 

.00037 

00001  00000 

7 

♦ 

-00015  - 00014 

00016 

00017 

< 000074 

000' 't, 3 

00032 

,00018  ,00015 

8 

,00002  00002 

,00001 

00000 

000173 

000172 

,00006 

„ 00007 

.00027  ,00026 

9 

„ 00007  ,00007 

00010 

,00010 

,000027 

,000031 

c 00001 

0 00001 

,00018  ,00018 

10 

,000025 

- 000023 

11 

0 000068 

.000063 

12 

,000012 

.,000013 

13 

.000011 

,000010 

14 

,000033 

,000034 

fc.O+oog  i“.05)39 

s 

Z*.  04/73 

/'j- 

'Z*.o+os/a 

^jD23oqc 

§_*.033TT 

/ 

^.o/szs 

/ 

* £ 

2 9.03/oj  2.S ,0/Z0! 
' / 

/ ' 

5*  ,04o/3  Sa‘°3f  i y 

S*o+)93 

S^.CZ&CiP 

$*.0Z3iS* 

.O/STU 

5-  a .f«  .0/233 
1 

«x  ■* 


101 


Table.  XI  v u;  ic.or.ti ) 
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Table  III 


The  Complex  Series  .*4'  * J Enter  i:«;  ;he 

First  Varii'.ticnal  Tern  for  2 , for  Various  4"? /?'. 
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Table  III 
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Table  IV < First  Variational  Approximation  to  Series  Impedance  for 

Change  of  Cross  Section  (!io  change  of  mate >p 

and  to  Reflection  Coefficient  /j?  (Semi -infinite  Termination) 
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T&bic  IV 


First  Variational  Approximation  to  Ses^iew  Impedance  for 
;-ha.ig<;  of  Cross  beet  ion  (No  change  of  material) 

and  to  Reflection  Coefficient  Q (Semi “inf inite  Termination) 

( & is  in  radians  and  Is  in  th<»  3rd  quadrant  i ,, e „ — 77-5 & ) 
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Table  V (.a)  Hither  Order  Series  d 
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Table  V (b)  'valuer  ol*  ■*/£,  and  & to  Second  Variational 

Approximation  at  tfte.je  selected  K and  jc  ‘fi* 1 
values  (vil.h  first  approximation  repeated) 


— 

(<%'/?  f 

r,m 

r 

<> 

l 

i 

1 

1 

/£/ 

G+  7T 

>0011 

t-Z 

A 00111  - 

-001111 

l 

1 

1796  | 

.0043 

B 00102  " 

0 001021 

.1796 

.0039 

Oil 

/ 2 

A 01109 

.011031 

,1723 

.0439 

B .01016 

.010191 

.1734 

-0404 

5o0i 

t.z 

A .03501  ~ 

,067291 
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